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. Abstract 

We discuss the ill-posed Cauchy problem for elliptic equations, which is pervasive in 
inverse boundary value problems modeled by elliptic equations. 

■ We provide essentially optimal stability results, in wide generality and under sub- 
stantially minimal assumptions. 

As a general scheme in our arguments, we show that all such stability results can be 
derived by the use of a single building brick, the three-spheres inequality. 
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1 Introduction 

> 

Hadamard, in his paper of 1902 [37] where he laid the basis of the notion of well-posed 
problems, used the Cauchy problem for Laplace's equation as his first example of a problem 
which is not well-posed. Later, in 1923, he published [38] his well-known example of instabil- 
ity, see Subsection II. II below for further discussion and also Maz'ya and Shaposhnikova [70] 

■ for additional information. 

It may be curious to note that in the same span of time, the physical relevance of this 
problem was to be encountered in the applications. For instance, in geophysical under- 
ground prospection, the geoelectrical method was initiated in those years, see for instance 
Stefanesco et al. [92] and the historical account in Zhdanov and Keller [100]. And in fact it 
is well acknowledged, by now, that the geoelectrical method involves, even in its most basic 
formulation, the solution of a Cauchy problem for Laplace's equation! On this respect, one 
can consult, for instance, the initial considerations in the book by Lavrent'ev, Romanov 
and Sisatskii [65]. 
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Nowadays, it is widely recognized that the Cauchy problem for Laplace's equation, 
and more generally for elliptic equations, has a central position in all inverse boundary 
value problems which are modeled by means of elliptic partial differential equations, Inverse 
Scattering, Electrical Impedance Tomography, Optical Tomography, just to mention a few. 
The continuing interest on this kind of problem is documented by the number of publications 
which are currently appearing on this problem. For instance, we have recorded at least 15 
papers explicitly devoted to this topic, which have appeared in the last three years on this 
Journal. 

Therefore, we believe that it might be useful to formulate in a clear fashion the state of 
the art on the issue of stability, which is obviously a crucial cornerstone of the convergence 
analysis of any reconstruction procedure and also of the stability analysis of many nonlinear 
inverse boundary value problems, whose treatment involves, in one way or another, the 
analysis of an ill-posed Cauchy problem. 

In this introduction we do not intend to present a complete discussion on the historical 
development on this subject since Hadamard, because various monographs, Lavrent'ev [64], 
Payne [83], Lavrent'ev, Romanov and Sisatskh [65], Isakov [47,48], Hormander [43,45], 
already contain abundant information on such development. However, besides stating and 
proving results of stability of sufficient generality and optimality it may be useful to discuss 
some different, although intertwined, lines of reasoning which, in our view, have led to the 
current state of the art. 

1.1 Instability and conditional stability 

In his essay of 1923 [38], Hadamard provided a fundamental example which shows that a 
solution of a Cauchy problem for Laplace's equation does not depend continuously upon 
the data. The example is as follows. 

Consider the solution u = u n , n = 1,2,... to the Cauchy problem in the upper half 
plane 

Au = 0, in {(x,y) G K 2 | y > 0}, 

u(x, 0) = 0, for every i£K, 

u y (x, 0) = A n sinnx, for every x 6 R. 

We have 

u n = — - sinnx sinhny. 
n 

If we choose A n = — or A n = — for some p > 0, or even A n = e _v/ ™, it turns out that 

n nP 

Un,y(x, 0) — ► uniformly as n — > oo 

whereas, for any y > 0, 

u n (x,y) = — sinnxsinhny blows up as n — > oo. 
n 

As is well-known since Tikhonov [95] the modern notion of stability for ill-posed problems 
(also called conditional stability) consists of estimating the dependence upon the data of 
the unknown solution of the problem at hand, when an a-priori bound on the solution itself 
is available. 
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It is a remarkable fact that Hadamard eventually acknowledged that continuous depen- 
dence can be restored in presence of an a-priori bound. In fact, in his treatise of 1964 [39, 
p. 146] he wrote 

"D'apres un remarquable resultat du a M. Pucci, l'absence de continuity de la 
solution u du probleme de Cauchy considered comme functionelle des donnees 
initiales est solidaire du fait que cette solution est susceptible d'augmenter in- 
finitement: M. Pucci constate que le choses changent si Ton connait une borne 
superiore de le valeur absolute |«|."B 

In fact, Hadamard is referring to a paper by Pucci [88] of 1955 where one of the first results 
of stability for the Cauchy problem for Laplace's equation was obtained. Let us recall that 
in the same years other stability estimates were obtained by John [52], Landis [60] and 
Lavrent'ev [62,63], see also Pucci [89] and John [53]. 

It may be instructing to observe that the same example by Hadamard may be used to 
exhibit the best possible rates of continuous dependence in presence of an a-priori bound. To 
this purpose let us describe a Cauchy problem in the most simple and favourable setting. We 
express the a-priori bound and the bounds on the data with respect to norms which can be 
considered as the natural ones in the standard variational formulation of Laplace's equation, 
but it will be evident that analogous results would be obtained also if other (reasonable) 
functional frameworks are considered. 

Consider a Cauchy problem in a rectangle 

{Au = 0, in (0,vr) x (0, 1), 

u(x,0) = 0, for every x £ (0,7r), 

u y (x,0) = ijj(x), for every x G (0,7r), 

and, in order to make things even simpler, we further assume a zero Dirichlet condition on 
the vertical sides of the rectangle 

(1.2) u(0, y) = u(-k, y) = 0, for every y G (0, 1). 

The natural function space of the solution u, in a variational setting, is the Sobolev space 
^((O.tt) x (0,1)) and thus, as a-priori information on the unknown solution u, we assume 
the bound on the Dirichlet integral 

(1.3) // (u 2 x + Uy) dxdy < E 2 , 

J J (0,7r)x(0,l) 

for a given E > 0. 

The prescribed inhomogeneous data if), which express the (partial) Neumann data on 
the lower horizontal side of the rectangle, naturally lives in the trace space H~ 2 (0, 7r). Let 
us assume then that the following error bound is known 

(1.4) \U\\ 1 < 77, 



* "After a remarkable result due to Mr. Pucci, the lack of continuity of the solution u of the Cauchy 
problem considered as a functional of initial data is joint to the fact that this solution is susceptible of 
infinite growth: Mr. Pucci observes that things change if one knows an upper bound of the absolute value 
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for some given rj > 0. 

A stability estimate would consist of a bound of some norm of u evaluated inside the 
rectangle in terms of some function u (rj, E) which should be infinitesimal as rj — > + . The 
Hadamard example provides us with limitations on such infinitesimal rate. Let us choose 
once more 

(1.5) ip n (x) = A n sin nx, n = 1, 2, . . . 

and let us select A n in such a way that equality holds in (jl.3p . We obtain 

(1-6) Al = ~-^-E 2 . 

7rsmh2n 

Consequently, in (|1.4[) we have equality when r\ = ry„ where rj n is given by 

(1.7) rj 2 = E 2 — 4E 2 e~ 2n , asn^co. 

sinh2n 

If we wish to estimate the L 2 -norm of u in the rectangle (0, tt) x (0, T), for some T £ (0, 1], 
then we see that the solution to (jl.ip . (|1.2p with ip given by (|1.5|) . (jl.6p . satisfies 

, . „ ll2 -E 2 /sinh2nT \ 2 e 2n( - T ~^ 

(1-8) lknlL 2(( o l7r )x(o,T)) = — ^ y 2n " 1 J ~ ^ ~2^-' as n ^ 00. 

That is 

(1-9) IMIl 2 ((o^)x(o,t))~^^J l, l0g ^J ,asn^oo. 

Therefore, if T < 1, then the stability of the determination of u up to the level y = T is 
at best of Holder type. Whereas, if we want to recover u in all of its domain of definition 
(up to the top side of the rectangle, where no boundary data is prescribed) then the best 
possible rate of stability is logarithmic. 

It is generally acknowledged that this phenomenon has a quite general character when 
dealing with the Cauchy problem for elliptic equations. And in fact the results in the 
following sections agree with such a scheme. 

We shall distinguish between two types of results. 

Stability estimates in the interior. A solution u of a Cauchy problem is a-priori known 
to be bounded (with respect to some norm) on a connected open set Q, Cauchy data 
are prescribed on some portion S of d£l and we wish to estimate u on some connected 
open subset G of Q which is at a positive distance from d£l \ S, the part of the 
boundary where no data are prescribed. 

Global stability estimates. We want to estimate u in some norm in all of f2 when an 
upper bound with a slightly stronger norm in the same set f2 is a-priori known, and, 
as before, Cauchy data are prescribed on a portion £ of d£l. 

We reiterate that the Hadamard example tells us that for a stability estimate in the 
interior we cannot expect anything better than a Holder rate, whereas for global stability 
the optimal rate will be of logarithmic type at most. 
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A rather general treatment of stability in the interior is due to Payne [82] (see also 
Payne [81], Trytten [96]). Global stability estimates in a wide generality, that is for general 
elliptic operators and general domains, have been known and used for quite a while, but 
probably statements and proofs are not easily available in the literature. In fact, a surge of 
interest on this topic occurred in the 90's in connection with nonlinear inverse boundary 
value problems with unknown boundaries and global stability estimates were described 
and used (Alessandrini [4], Alessandrini and Di Benedetto [6], Beretta and Vessella [16], 
Bukhgeim, Cheng and Yamamoto [20,21], Alessandrini, Beretta, Rosset and Vessella [5], 
Cheng, Hon and Yamamoto [26,27]) but, unfortunately, most of the times such estimates 
were not explicitly stated as independent results. As an exception, we mention Takeuchi 
and Yamamoto [94, Theorem 10]. 

1.2 Analytic continuation 

In the special case of two space variables, the Cauchy problem for Laplace's equation is 
equivalent to the problem of continuation of a complex analytic function from values pre- 
scribed on an arc. This problem was treated with great ingenuity by Carleman in 1926. His 
theory is expounded in the essay of 1926 [24] and his ideas have had a great influence in the 
subsequent developments of the theory. A modern treatment of the connection between the 
Cauchy problem and the analytic continuation in two dimensions can be found in [65, Ch. 1, 
§ 2]. The crucial tool for stability is the so-called method of harmonic measure. In fact, the 
seminal idea of this method can be traced back to Carleman [23], but a general formula- 
tion of this approach can be attributed to F. and R. Nevanlinna [79], see in this respect 
Goluzin [35, Ch. VII, § 4]. 

The two-dimensional theory for Cauchy problems maintains a special position also when 
dealing with elliptic operators with variable coefficients. Indeed, still with the aid of complex 
analytic methods, uniqueness, Alessandrini and Magnanini [8], and stability, Alessandrini 
and Rondi [10], can be obtained with no need of regularity assumptions on the coefficients 
(contrary to what happens in higher dimensions, see in this respect Subsection 11.31 below) . 
An extended version of the harmonic measure technique has been developed also for the 
variable coefficients case in Alessandrini and Rondi [10]. In Section [3] below we shall discuss 
in more detail this kind of results. 

1.3 The Cauchy problem and the unique continuation property 

An issue which is strictly related to the Cauchy problem is the one of unique continuation. 
An elliptic operator C is said to have the (weak) unique continuation property if for any 
solution u to Cu = in a connected open set O C M n , that vanishes on an open subset 
G C O, it follows u = in S7. A general proof of the equivalence of the uniqueness of 
the Cauchy problem with the weak unique continuation property can be traced back to 
Nirenberg [80]. This equivalence has been especially important in establishing the limits of 
validity of uniqueness in terms of the regularity of the coefficients of the elliptic operator 
involved. In fact it was shown by Plis [86] that for an elliptic operator C in dimension n > 3, 
the unique continuation may fail if the coefficients of the principal part are Holder continuous 
of any exponent smaller than 1. Further examples were obtained by Miller [71], see also 
Miller [72] for improvements and discussion on such counterexamples, and Mandache [69], 
Filonov [31] for further developments. 
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The progress on this issue of unique continuation was initiated by Carleman [25], subse- 
quent advances were due to Miiller [76] , Heinz [42] , Hartman and Wintner [40] , Cordes [28] , 
Aronszajn [14]. Eventually, it was proved by Aronszajn, Krzywicki and Szarski [15] that 
the unique continuation property holds true when the coefficients in the principal part are 
Lipschitz continuous. Soon afterwards, Plis produced the already mentioned example that 
shows that such Lipschitz continuity provides indeed the crucial threshold. A great deal 
of investigation followed, especially with the purpose of extending the unique continuation 
property to Schrddinger operators with singular potentials. With no ambition of complete- 
ness, let us mention Hormander [44], Jerison and Kenig [50], Garofalo and Lin [32,33], 
Fabes, Garofalo and Lin [30] and, more recently, Koch and Tataru [55]. 

1.4 Three-spheres inequalities 

The unique continuation property is also connected to the problem of the stability for the 
Cauchy problem. In fact, proofs of the unique continuation property depend on inequalities 
which can also be applied to the estimation of stability. There are two families of such 
inequalities 

• Carleman estimates, 

• Three-spheres inequalities. 

Both types have been succesfully used in the study of stability, and they are strictly 
intertwined, in fact three-spheres inequalities can be deduced by Carleman estimates. 

In this paper, a central theme that we intend to stress is that three-spheres inequalities 
can be used as a universal building brick to derive optimal stability estimates. Our guiding 
idea shall be to confine all the hard-analysis which is required to the derivation of a basic 
inequality in a simple geometrical setting (the three-spheres) and then use it iteratively to 
adapt to general geometrical configurations. 

We find quite instructing to remark at this point that Hadamard, the same mathemati- 
cian who first pointed out the ill character of the Cauchy problem, was the one who first 
provided a cure (actually, even before the illness was diagnosed!). In fact, Hadamard first 
stated in 1896 [36] a three-circles inequality, which in its simplest manifestation is as follows. 
Given a holomorphic function / in the disk {z E C | \z\ < R} then the function 




is convex. Quoting once more Hadamard [36, p. 94] 

". . . je designerai par r\ le logarithme du module maximum de la fonction sur le 
cercle de rayon (oil £ est un nombre reel quelconque). Le lieu du point (C,v) 
est une courbe C qui tourne toujours sa concavite vers le r\ positifs; ... "0 

' "... I shall denote with r\ the logarithm of the maximum modulus of the function on the circle of radius 
(where f is any real number). The locus of the points (£,r?) is a curve C that always bends its concavity 
towards the positive rfs; ..." 
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In other terms, if < r\ < r2 < r^ < R, then 

max |/ (z)\ < (max \f (z)\) (max \ f (z)\) 

\z\=r 2 \\z\=ri / \N=»*3 / 

where a £ (0, 1) is given by 




This inequality had a great influence in the following development of complex analysis [29] 
but it had also a seminal character in the study of unique continuation for elliptic equations. 
A three-spheres inequality for elliptic operators whose principal part coefficients are C 2 was 
proved by Landis [61], and in fact his proof was based on Carleman's type estimates. Another 
proof obtained by a method of differential inequalities for integral norms, which took the 
name of logarithmic convexity, was obtained by Agmon [3]. In the 70's the general concept 
of logarithmic convexity had indeed a notable influence in the analysis of various ill-posed 
problems for partial differential equations. The proceedings book edited by Knops [54] 
documents the advances in this direction. 

With more precise connection with the stability for the Cauchy problem, this approach 
reached its apex in the work by Payne [81,82], see also the contribution by Trytten [96]. 
In particular in [82], Payne obtained stability estimates in the interior of Holder type when 
the coefficients of the principal part are C 1 . 

Returning however to the more specific issue of the three-spheres inequalities, we men- 
tion that recent proofs have been obtained by Brummelhuis [19], and by Kukavica [58]. 
Both authors use, with some variations, the so-called method of the frequency function by 
Garofalo and Lin [32] . In fact the frequency function method can be viewed as a further 
notable advance and clarification of the ideas of logarithmic convexity. 

We shall formulate various versions of the three-spheres inequality and we shall provide 
a proof which is modeled, with few adjustements, on the one by Kukavica [58]. Let us 
remark however that the same inequality might be obtained also by means of Carleman 
estimates. We refer to Vessella [99] for a general discussion of this approach, in the wider 
context of parabolic equations. Let us also quote Lin, Nakamura and Wang [67] for a very 
recent investigation in this direction. 

1.5 Doubling inequalities 

It cannot be omitted at this point that reasearch on unique continuation has been especially 
concentrated on the aspect of the strong unique continuation property, that is, if u solves 
the elliptic equation Cu = in a connected open set and u vanishes of infinite order at 
one point xo G SI (that is u (x) = O (x — xo\ N ^j as n xo, for every N = 1, 2, . . .) then u 
has to be zero everywhere. 

This property does not have a direct connection with the stability of the Cauchy problem. 
However, quantitative versions of the strong unique continuation property have shown to 
be very useful in the study of stability of certain inverse boundary value problems. Such 
quantitative estimates are in fact the doubling inequality by Garofalo and Lin [32] and the 
doubling inequality at the boundary, see Adolfsson, Escauriaza and Kenig [2], Adolfsson 
and Escauriaza [1], Kukavica and Nystrom [59]. Applications to elliptic inverse boundary 
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value problems occurred for instance in Alessandrini, Rosset and Seo [12], Alessandrini, 
Beretta, Rosset and Vessella [5]. We recall also that such quantitative estimates have been 
used in connection with the problem of estimates of continuation from measurable sets, 
which is very much related to the Cauchy problem as well. In this respect let us mention 
Nadirashvili [77,78], Vessella [97,98], and Malinnikova [68]. 

1.6 The scheme of a stability proof 

Our general scheme of proof for the stability of a Cauchy problem will be as follows. 

I) First we prove a three-spheres inequality. We shall need an inequality of this kind not 
only for solutions of homogeneous elliptic equations but more generally for solutions of 
inhomogeneous equations Cu = T with an H^ 1 right-hand side J- (see Theorem ll.lOp . 

II) Next we use iteratively the three-spheres inequality to obtain estimates of propagation 
of smallness. We assume an a-priori bound on a solution u on its domain f2 and that 
u is small in a given (small) ball B ro (xq) C O and we estimate how small is u in some 
larger connected open set G C ft. If G is at a positive distance from d£l we shall speak 
of estimates of propagation of smallness in the interior (see Theorem 15. ip , instead if 
G agrees with f2 we shall speak of global estimates of propagation of smallness (see 
Theorem 15.31 and also Remark 17. 2p . 

III) Then, given a solution of a Cauchy problem in Q with data on a portion E of <90, we 
extend u to an open set A outside of f2, whose boundary agrees with dfl on a subset of 
E. We perform such extension in such a way that the extended function u solves in A 
an inhomogeneous equation with a right-hand side which is controlled in a Lipschitz 
fashion by the Cauchy data on E (see Theorem 16.21 and the following Remark 16. 3p . 

IV) Finally, we apply the estimates of propagation of smallness in the augmented domain 

o 

Q = f2 U A thus obtaining an interior stability estimate for the Cauchy problem in 
Theorem 1 1.7\ and a global stability estimate in Theorem 11.91 

1.7 The main hypotheses and statements 

In order to discuss in more detail the hypotheses that are used in our treatment it is 
necessary to introduce some notation and definitions. 

Given x G R n , we shall denote x = (x ! , x n ), where x' = (xi, . . . , x n _i ) G R"- 1 , x n G K. 
Given x G R n , r > 0, we shall use the following notation for balls and cylinders. 

B r (x) = {y G R n | \y-x\< r}, B r = B r (0), 

B' r (x') = y G R"- 1 I \y - x'\ < r}, B' r = B' r (0), 

T a ,b(x) ={y = (y 1 , y n ) € R n | \y' - x'\ < a, \y n - x n \ < b}, r ajb = 1^(0). 

Throughout this paper we shall denote by il a bounded open connected subset of R n . 
In places we shall assume that the boundary of $7 is Lipschitz according to the following 
definition. 
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Definition 1.1 (Lipschitz regularity) We say that the boundary of f2 is of Lipschitz 
class with constants po, Mq > 0, if, for any P E d£l, there exists a rigid transformation of 
coordinates under which P = and 

(1.10) nnT^ p AP) = {x = (x',x n )£Tp iL I x n >Z(x')}, 
where Z : B'^_ —* R is a Lipschitz function satisfying 

A/Q 

(1.11) Z(0) = 0, 



(1.12) ||Z|| / x < M oPo . 



Win' 

V W 

Remark 1.2 For practical purposes it will turn out useful to assume throughout that 
Mo > 1. In fact it is evident from Definition 11.11 that conditions (|1.10p — (|1.12|) continue to 
hold if Mo is increased. 

Remark 1.3 Throughout this paper we shall use the convention to normalize all norms in 
such a way that they are dimensionally equivalent to their argument and coincide with the 
standard definition when the dimensional parameter po equals 1. For instance, the norm 
appearing above is meant as follows 

\\ z \\ ( \= ll-Zll / \ +Po||VZ|| 

C°.l B> L°° B' L°° I B' 



Similarly, we shall set 



Mli 2 (n) = Po 2 



u 2 



i 

2,2/ IV7..I2 1 ~ 



\\u\\m(u.) = Po 2 u + po / |Vu 
\Jn Jn 

and so on for boundary and trace norms such as II • II r 2 r«ox , II • II 1 , II • II 1 

In some instances we shall require that only a limited open portion E of 517 be Lipschitz 
in the following sense. Some further notation is necessary We shall denote 

(1.13) E' = 0n\E, 
and, for every P G E, we set 

(1.14) r{P) = dist(P,E'), 



(1-15) p(P)=minL , r(P)M ° 



vr+Mj 
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Definition 1.4 An open subset E C 80, is said to be an open Lipschitz portion of 90 with 
constants po, Mq > 0, if, for any P G E, there exists a rigid transformation of coordinates 
under which P = and 

(1.16) finr^ #) (P) = {x = (x',x n ) er^ #) | x n > z(x')}, 

where Z : P/ P0 — > K is a Lipschitz function satisfying (|l.lip - (|1.12p . 

We shall also need an assumption on E in order to control jrom below its smallness. 

Definition 1.5 We shall say that S has size at ieasi pi, < p\ < po, if there exists at least 
one point P G E such that 

(1.17) p{P)>Pi- 

The elliptic operators that we shall consider are of the following form 

(1.18) Cu = div (AVu) + cu 

where A = A(x) = {a,ij(x)}, x G R n , is a real- valued symmetric n x n matrix such that its 
entries are measurable and it satisfies, for a given constant K > 1, the ellipticity condition 

(1.19) K' 1 ^ 2 < A(x)i ■ f < K\i\ 2 , for almost every x G R n , for every ^ G M n . 

Furthermore, when n > 3, we also assume that, for a given constant L > 0, the following 
Lipschitz continuity holds 

(1.20) - A(y)\ < — \x - y\, for every x,y G W 1 . 

Po 

Concerning the zero order term, we assume that c G L°°(M. n ) with 

/ \ 1 1 1 1 ^ 

(1-21) ||c||,L°o( R n) < — • 

Po 

Here and for the rest of this paper we shall assume that A is a symmetric matrix 
of coefficients satisfying the ellipticity condition (|1.19p and such that, if n > 3, it also 
satisfies the Lipschitz condition (ll.20p . We emphasize that, in all the following statements, 
whenever a constant is said to depend on L (among other quantities) it is understood that 
such dependence occurs only when n > 3. 

We remark once and for all that also first order terms with bounded coefficients could be 
added to the treatment with little additional effort. We have chosen to confine ourselves to 
the above variational structure (|1.18p because we believe that (|1.18p provides a sufficiently 
wide and useful setting for enough applications. We insist however that we allow the presence 
of the zero order term, with no sign, nor smallness, condition on the coefficient c because 
of the importance of the applications to wave phenomena at a fixed wavenumber. 

Let us now introduce the rigorous weak formulation of the Cauchy problem that we 
shall use. First it is necessary to introduce some further notation and some function spaces. 

The solution u will be assumed to belong to the H 1 (Q) space and the following a-priori 
bound will be prescribed 

(1.22) \H\ mm <E, 
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for some given E > 0, for the purpose of a global stability estimate, whereas for the stability 
in the interior we shall more simply require that, given Eq > 0, 

(1.23) \\u\\ L2(n) < E . 

We shall fix E as an open connected portion of d£l and we shall always assume it to be 
Lipschitz. 

We shall consider as test functions space the space -ff^ftUE) consisting of the functions 

I 

ip G H 1 ^) having support compactly contained in OUS. We denote by -£f c 2 (E) the class 

of i?2 (E) traces of functions tp G Hl (Q U E). We then define H~z (E) as the dual space to 
1 

-ff c 2 o(£) based on the £ 2 (E) dual pairing. 

The Cauchy data g, ip will be taken in their natural trace spaces, namely 

(1.24) fl eff5(E)^ £ i[4(E). 
As bounds on the Cauchy data we require 

(1-25) M H i {s) +PoUW H - ks) <V- 

To begin with, let us consider the following more or less standard formulation of a 
Cauchy problem 



(1.26) 



div (AVu) + cu = f, in ft, 
u = g, on E, 

AVu ■ v = on E, 



where the right-hand side / can be assumed in 1? (W 1 ). 

The corresponding rigorous weak formulation of the Cauchy problem would be to find 
u £ H 1 (f2) such that U| S = g in the trace sense and 

(1.27) / (AVu ■ Vip - cucp) = ipip- ftp, for every ip £ H x co (O U S) . 

Here, the integral J s ipip is to be properly interpreted as the dual pairing between H~2 (X) 

i 

and Hco (X). Note that the right hand side in (jl.27j) represents a bounded linear functional 

over ^ Q (QUS). An even more general bounded functional over H^QUY,) could be written 

as follows 

(1.28) Hv) = [ ipv- [ (f<p-F- Vip) , for every p G R\ (ft U E) 

where F G L 2 (M n ;l n ) is a given vector field with 1? components. It will be convenient for 
us to admit such type of right-hand side in our formulation. Observe, however, that such a 
representation of members of (H^ Q (Q U E))* is exceedingly non- unique (and indeed we shall 
fruitfully take advantage of this non-uniqueness in Theorem I6.2p . 
Hence we shall prescribe 

(1-29) ll/llz,2(K«) + — l|-F1ll,2(R™ ; Rr*) < — , 

F0 Pq 

for a given e > and we shall formulate the rigorous weak formulation of the Cauchy 
problem as follows. 



11 



Problem 1.6 (The weak formulation of the Cauchy problem) To find u G H l (Q) 
such that = g in the trace sense and 

(1.30) / (AVu ■ Vip - cuip) = ipp- (ftp — F ■ Vtp) , for every <p G U E). 

We remark that an abstract interpretation of the above stated problem could be expressed 
as follows 



(1.31) 



div (AVu) + cu = T, in (#* (^UE)) S 
u\z=g, miJ5(S), 



where .F is given by (jl.28p . Unfortunately, this abstract formulation hides within its first 
equation a boundary contribution which should express the Neumann condition on E. 
Roughly speaking, such boundary contribution can be detected, if one formally integrates 
by parts the terms involving V<p appearing in (jl.30p . and formally obtains 



div (AVu) + cu = f + divF, in O, 
u = g, on E, 

AVu ■ v = ip + F ■ v, on E, 



and it is evident that this formal expression incorporates as a special case (|1.26p . We re- 
iterate, however, that the true interpretation of the Cauchy problem that we shall use is 
given by Problem 11.6^ but at the same time we emphasize that the full strength of such a 
rigorous formulation shall be used at a single step in our arguments only, and specifically 
in the above mentioned Theorem 16.21 In all remaining estimates occurring in this paper we 
shall merely make use of the interior weak formulation of the elliptic equation 

div (AVu) + cu = f + divF, 

which can be viewed in the more customary sense of i^ _1 (0) = (Hq(Q))* . 
Our main stability estimates are contained in the following theorems. 

Theorem 1.7 (Stability in the interior for the Cauchy problem) Letu G H 1 ^) be 
a weak solution to the Cauchy Problem \1.6l where E satisfies the conditions in Defini- 
tion\L]\and Definition\Tjh f G L 2 (R n ) and F G L 2 (R n ;M. n ) satisfy (fi~29l) and g G (E), 
if) G H~z (E) satisfy (ll.25p . There exists h, < h < with only depending on K, L, 
k, Mq and such that, assuming the a-priori bound 

(1-32) HL>(n)<^, 

then, for every h, < h < h, and for every connected open set G C SI such that 

(1.33) dist(G, dfl) > h, 

(1.34) dist(P,G)< l>] 



8M 
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where P G £ is the point appearing in Definition \1.5l we have 



(1.35) \\u\\ L 2 {G) < C[e + V ) S (E + e + r?) 1 ^, 

where C > and 5 E (0,1) satisfy 



i 



(1.36) C = C ^§ 

and 

Coin 

(1.37) <5 > a~^~ 

with a G (0, 1) only depending on K, L and k, C\ only depending on K , L, n, Mq and 
and C2 only depending on K . 

Remark 1.8 Let us observe that stability with a Holder rate could be obtained also when 
the connected subset G of Vt is allowed to touch the boundary portion E, while remaining 
at a positive distance from its complement X'. This fact will turn out to be evident from 
the proof, see also Remark 16.51 We have chosen the present formulation because in this way 
a more effective evaluation of the constants C and 8 is obtained in terms of the parameter 
h, which controls the distance of G from the whole boundary dQ. It would also be possible 
to obtain analogous results in terms of the distance of G from £' = d£l \ X at the price of 
assuming more information on the shape and regularity of the boundary of S within dil. 

Theorem 1.9 (Global stability for the Cauchy problem) Let Q be a connected open 
set of Lipschitz class according to Definition li.il Let u £ H 1 ($7) be a weak solution to the 
Cauchy ProblemEJb where f G L 2 (R") and F G L 2 (M. n ;M n ) satisfy ([09]) and g G (E) ; 
ip G H~2 (£) satisfy (ll.25|) . If u satisfies the a-priori assumption 

(1-38) IMlHipi) < E, 

then 

(1-39) IMI^n) <{E + s + r))u [ E £ + e \ J , 



where 

C_ 



(1.40) u{t)<- T - JI , fort<\. 



where C > and fi, < fj, < 1, only depend on K, L, k, Mq, ^ and 

According to the scheme already illustrated in Subsection 11.61 our first step will be the 
proof of a three-spheres inequality. In fact our basic building brick will be the following 

Theorem 1.10 (Three-spheres inequality) Let u G H 1 (Br) be a weak solution to the 
inhomogeneous elliptic equation 

(1.41) div(AVu) +cu = f + divF, in B R , 
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where f G L 2 (R n ) and F G L 2 (R n ;~R n ) satisfy (jl.29p . Then, for every t\, r2, r^, with 
< rx < r 2 < r 3 < R, 

( L42 ) \M\L*(B r2 ) < C (\\u\\ L 2 {Bri) + e) (||u||L2(B r3 ) +ej 

where C > and a, < a < 1, on/y depend on K, L, k, max|j|, l\, ^ and 

Our strategy for the proof of such an inequality breaks down in the following steps. 

a) First we prove a three-spheres inequality for the homogeneous equation in pure principal 
part (Theorem 12. 11) with some limitations on the radii, namely < n < r 2 < r^/K < 
r 3 < R, when n > 3. In two dimensions, a three-spheres (in fact circles) inequality may 
be obtained, with a different technique, without the limitations on the radii and with a 
possibly discontinuous coefficient matrix A (Theorem 13. lip . 

b) Next we adapt the proof to operators also containing the zero order term (Theorem l4.ip . 

c) We include the presence of the right-hand side in the equation (Theorem 14. 4p . 

d) We remove the limitations on the radii, finally obtaining the above stated Theorem ll.101 

We shall also consider estimates of propagation of smallness, in Section The main 
results will be Theorem 15.11 (propagation in the interior) and Theorem 15.31 (global propa- 
gation). In fact, Theorem 15.11 is obtained as a consequence of the three-spheres inequality 
with restrictions on the radii, Theorem 14.41 and the above mentioned step|d]) is obtained as 
a special case of Theorem 15.11 

In Section [6l we prove the extension argument (Theorem I6.2p which enables to apply 
the estimates of propagation of smallness (Theorems 15.11 and 15. 3p to complete the proofs of 
Theorems 11.71 and 11.91 

Finally, we shall address the problem of global stability when no regularity assumption 
is available on the domain O. This situation in fact often occurs in inverse problems with 
unknown boundaries [4, 5, 10, 22, 73, 74, 91, 99]. If two solutions of the same equation have 
different domains of definition, then their difference solves an equation in the intersection 
of the two domains. Such an intersection may be, in principle, highly nonsmooth, even 
if the starting domains satisfy some a-priori regularity assumptions. It is then useful to 
derive a preliminary, maybe very weak, global stability for such difference of solutions in 
their common domain of definition (or at least in one of its connected components). We 
believe that the present formalization of this argument, expressed in Theorem 17. 11 might be 
a useful tool for future use in other inverse problems with unknown boundaries. Of course 
some assumptions will be needed anyhow. The basic one that we shall use on Q is rather 
weak and can be summarized as follows. We shall assume that there exists a family {Gh} of 
connected subsets invading f2 such that the measure of the difference | f2 \ Gh | is controlled 
by a given power h of the distance h of Gh from <9f2. It may be easy to predict that, with 
such weak hypotheses, only a weak result can be obtained. In fact we are able to achieve 
only a stability of log-log-type. 
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1.8 Concluding remarks 

We believe that the above theorems provide optimal results of stability for the Cauchy 
problem. The optimality is achieved in many respects. 

1) The regularity assumptions on the coefficients of the principal part are minimal, Lipschitz 
continuity when n > 3 and L°° when n = 2. In fact these are known to be optimal 
conditions for uniqueness. 

2) The regularity assumptions on the boundary are kept to a minimum. In fact assuming 
Lipschitz regularity of the portion £ of dQ where the Cauchy data are assigned seems 
to be a nearly minimal condition, just in order to give sense to the Cauchy problem. On 
the other hand, when dealing with global stability, some conditions on the boundary, 
which may ensure the uniform reachability of the boundary points from the interior, 
seem necessary and Lipschitz regularity appears to be enough general and at the same 
time meaningful in terms of applications. More general reachability conditions could 
be considered, such as the NTA corkscrew condition of Jerison and Kenig [51]. Let us 
mention in this direction the related investigation by Rondi [91]. 

3) The Cauchy data are evaluated in their natural spaces. This is indeed a slight improve- 
ment with respect to previous studies where typically i/ 1 (S)-norm on g and £ 2 (S) -norm 
on ifi are considered, see Payne [82], Isakov [48], Takeuchi and Yamamoto [94]. 

4) As already discussed, the stability rates obtained have an optimal character. The quanti- 
ties in the stability estimates which might require further investigation are the constants 
and the exponents (for instance C and 5 in Theorem 1 1 . 71 and C and fi in Theorem ll.9p . It 
would be interesting, although possibly rather difficult, to simultaneously optimize such 
pairs of quantities with respect to the geometry and also with respect to the coefficients 
of the equation. In connection to this issue we wish to mention some remarkable results 
by Hrycak and Isakov [46] and Subbarayappa and Isakov [93] who have considered the 
Helmholtz equation Au + k 2 u = and proved that the stability estimates improve as 
\k\ increases, in the sense that, as \k\ — > oo, the logarithmic term becomes negligible, 
whereas a Holder term prevails. See also Isakov [49] for a related analysis with a more 
general equation. 

The plan of the paper is as follows. Section [2] contains the proof of a three-spheres 
inequality for a homogeneous equation in pure principal part, the main result being Theo- 
rem [27TJ In Section [3] we investigate the two-dimensional case. We illustrate the connection 
with quasiconformal mappings in Proposition 13.11 We introduce the notion of /2^-harmonic 
measure, Definition 13.51 an d we apply it to stability estimates in the interior for Cauchy 
problems, Theorems 13.71 and 13.81 Finally, we deduce various forms of three-spheres (circles) 
inequalities, Theorem 13.101 and Theorem 13. 1 II which will be the one we shall use in the rest 
of the paper (together with Theorem 12. 1 j) . In Section 0] we adapt the previously obtained 
three-spheres inequalities in order to encompass equations containing the zero order term 
and an inhomogeneous right-hand side. The final result of this section is given in Theo- 
rem 14.41 Section [5] is devoted to estimates of propagation of smallness. In Theorem 15.11 
we prove the interior estimate. As we already mentioned, as a corollary we also derive the 
proof of the above stated Theorem ll.lOl Next, in Theorem 15 .31 we prove the global estimate. 
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Section [6] contains the proof of the two main stability theorems for the Cauchy problem, 
Theorems 11.71 and 11.91 As we already illustrated, Section [7] deals with a global stability 
estimate for the Cauchy problem when no regularity is assumed on Q, Theorem 17.11 



2 The basic three-spheres inequality, arbitrary dimension 

Let us begin by considering the homogeneous elliptic equation in pure principal part 

(2.1) div(AVu) = 0, in B R , 

where A satisfies (|1.19p and (jl.20p . for given constants K > 1 and L > 0. 

Theorem 2.1 (Three-spheres inequality — pure principal part) If the above stated 
hypotheses hold, for every t\, Ti, with < r\ < r2 < § < r$ < R, 



(2-2) \\u\\ L 2 {Br2 ) < Q\H 



LHB ri )\\ u \\ L 2m) 



\l-a 



where Q>1 only depends on K, L, max|j|,l|, and where 

(2.3) a = 

log^ + Clog^f' 

with C > only depending on K , L and max l|. 

Remark 2.2 It should be stressed here that, in the following sections, we shall not use the 
full strength of Theorem 12.11 with the essentially explicit formula (|2.3p for the exponent of 
the three-spheres inequality. Actually, we shall merely use the fact that a only depends on 
K, L, max j^,l j, j^-We have just recorded here the representation (|2.3|) because it 
may be of independent interest. We recall, in particular, that a three-spheres inequality (|2.2p 
with exponent in the form (|2,3p actually implies the strong unique continuation property 
in a rather straightforward fashion. The argument is as follows. We fix r2,r% and we allow 
r\ — * 0, in this case 

a = a(n) > - — X - u -. 

With no loss of generality we may assume also ||w||^2(^ r3 ) < 1. If for every N = 1,2,... we 
have |M|£2( Br ) = 0(r ) as r — > 0, then, for every N, there exists r' > such that 

\\u\\ L 2 {Bri) < (-^-j , for every r x < r' , 

consequently, by (|2.2|) . 

\\u\\ L 2 (Br2) < Q f^y° g7 ^ = Qe~ ClN - 0, as N - oo. 
See also Morassi, Rosset and Vessella [75] where this aspect was investigated in more depth. 
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As a first step, we shall prove a three-spheres inequality under the following additional 
normalization hypothesis 

(2.4) A(0) = Id. 

Given a weak solution u E H 1 {Br) to equation (|2.ip . let us define 

yl(x)x • x 



(2.5) /z(x) 



Ixl 2 



2 



(2.6) H{r) = / 

J 8 B r 

Theorem 2.3 (Three-spheres inequality — normalized case) If the previously stated 
hypotheses hold, for every r\, r2, r^, with < r\ < r2 < r$ < R, 

(2.7) H(r 2 ) < QH(r 1 ) tt H(r 3 ) 1 - a , 
where Q>1 only depends on K, L, max|^,ll, and where 



(2.8) a 



V2 

log a + clog * : 



with C > only depending on K , L and max 1 j. 
Remark 2.4 Since i^ -1 < n < K, we have 



(2.9) / . 2 <Q( / . 2 ) ( / u 2 ) 

J9B r2 \JdB ri J \JdB r3 J 

and, by integration and by Holder inequality, we also obtain 

(2.10) f u 2 <q( I u 2 \ ( j J\ . 

J Br 2 \J B r± J \yJ Br^ J 

Let us first derive Theorem 12.11 from Theorem 12.31 
Proof of Theorem 12.11 Let us introduce the change of variables 

(2.11) y = Jx, 

where J = \J .A _1 (0) and let us consider, for any r > 0, the ellipsoid 

(2.12) E r = {x e R n | (^(0) _1 x) • x < r 2 }. 
We have that £ r = J~ 1 (B r ) and 

(2.13) c Er C B^ r . 
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The function v(y) = u(J 1 y) satisfies the elliptic equation 

(2.14) div(AVu) = 0, in Br, 

Vk 

where A(y) = JA(J~ 1 y)J. It is straightforward to verify that 

(2.15) K- 2 \£\ 2 < A{y)£ ■ i < K 2 \i\\ for every y G M. n , £ G R n , 

(2.16) \A( yi ) - A(y 2 )\ < Kl — \yi - y 2 \, for every yi ,y 2 G E n , 

Po 

(2.17) 1(0) = Id. 

Therefore A satisfies the hypotheses of Theorem 12.31 Let p\ = -^=, p 2 = VKr 2 , p^ = -^|=. 
Since < n < r 2 < ^ < r% < R, we have that < p\ < p 2 < p3 < -7=, and by (|2.10p we 
have 



a / \ 1—a 

2 ^ r> I I „.2 1 ( / „.2 I 



(2-18) / v 2 <Q\ / v 2 \ v 

Jb P2 \Jb Pi J \Jb P3 J 

where Q > 1 only depends on K, L, max 1 1, and where 

log 

(2.19) a 



P2 



log £& + Clog ^' 

to P2 b Pi 



with C > only depending on i^T, L and max j^, 1 j. Coming back to the old variables 



we 



have 



(2.20) f u 2 <Q \ [ u 2 ) If 



a / \ 1—a 

2 I 



1/ 

P2 / V t, />3 



and recalling (|2.13p . the thesis follows. □ 

The proof of Theorem 12.31 is essentially based on the following three lemmas. 

Lemma 2.5 The function H(r) is absolutely continuous in (0, R) and, for almost every 
r G (0, R), 

(2.21) H'(r) = -[ {-Av ■ v + tr(A) +rd Xi a ij v j )u 2 + 2 ( uAVu-v. 

r JdB r JdB r 

Proof We have that 

^ 1 r 

u 2 A(x)x ■ v = 



(2.22) H(r) = - [ u 2 A{x)x ~ = - / 
r JdBr \A r J d 



•'1 r ^SB r 

- / div (tt 2 A(a;)a;) = - / ds / div (ti 2 A(x)x), 
r J B r r Jo JdB a 
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where v denotes the outer unit normal to B r . Since, by well-known regularity results [34, 
Theorem 8.32], u E C]^{B R ) and since A £ C^^Br), we have that div(u 2 A(x)x) £ 
L°°(Br) and therefore H(r) is absolutely continuous in (0, R). From (|2.22p . we have 



(2.23) H'(r) = -~H(r) + - f dW(u 2 Ax) = 
r r J dBr 

= -- / [in 2 + - / (tr(A))u 2 + - / u 2 d Xi a,ijXj + - / uu Xi a,ijXj, 
r JdB r r J dBr K r J dBr r J 9Br 

and, recalling that A is symmetric, (|2.21|) follows. □ 
Let us set 



(2.24) I(r) = / uAVu v= AVu ■ Vu > 0. 

Lemma 2.6 There exists a positive constant C, only depending on K and L, such that, for 
almost every r £ (0, R), 



(2.25) 



71 — 1 

H'( r ) - - — -H{r) - 2J(r) 



< —H(r) 



(2.26) I'(r)>2 [ -(AVu-v) 2 + ^^I{r) - —I{r), 

JdBr V r po 

Proof For almost every r G (0,-R), we can compute 

n — 1 

(2.27) H'ir) - - — -H(r) - 21 (r) = 

r 

= - I (—Au ■ v + tr(A) + rd Xi aijVj)u 2 — — / fiv 2 = 

r JdBr r JdBr 

= - I (-A(O)u ■ v + tr(A(0)))u 2 - — - / pu 2 + - / rd Xi a ij u j u 2 + 

r JdBr T JdBr r JdBr 

+ - f [-(A - A(0))v ■ v + (tr(A) - tr(A(0)))]u 2 = 

r JdBr 

n-l f s o , 1 f o 2 



(1 - n)u + - rd Xi aijVjU + 

r -IdBr r JdBr 



+ - f [-(A- A(0))u ■ u+ (tr(A) - tr(A(0)))]u 2 

'dBr 



r 



By (fl~20]) and we have that 



(2.28) |(A(x)-il(0))|<-|x| 



(2.29) II - u(x)\ = \(A(0) - A(x))v ■ v\ < —\x\ 

Po 



19 



and that the absolute values of the three terms in the right-hand side of (|2.27p are bounded 
by CH(r), obtaining (jZHj) . 
Let us notice that 



(2.30) I(r) = [ ds [ AVu • V 

Jo JdB a 

is continuously differentiable and that 



(2.31) I'(r) = / iVn ■ Vn. 

JdB r 

Let us recall the following generalization of the Rellich identity [90], due to Payne and 
Weinberger [84] 

(2.32) / aijd Xi ud X] uf k v k = 
JdBr 

= 2 fia k jd Xt ud Xj uv k - 2 fkO Xk ud Xi (aijd Xj u)+ 

J 0B r J B r 

+ / (d Xk f k aij - 2d Xk fia kj + f k d Xk aij)d Xi ud x 

J Br 



Br 

which holds for every / G C 0,1 (B r ;M. n ) and every u E H 2 (B r ), where aij are as above. 
Choosing 

< 2 - 33 » /w - ft- 

from (I238|) and (12391) . it follows that 

C 



(2.34) 



< 
Po 



where C is a constant only depending on L. Inserting (|2.33p in (|2.32|) and recalling (|2.3ip 
and (|ZMD , (|236]) easily follows. □ 

The proof of Theorem 12.31 shall be based on a differential inequality for the so-called 
frequency function, a notion first introduced by Almgren [13], see also [32]. 

(2.35) N{r) rUr) 



H(rV 



which is well-defined provided H{r) > 0. Observe that the thesis of Theorem 12.31 is trivial if 
u is identically constant in Br, or, as is the same, if I(R) = 0. Thus we may assume, with 
no loss of generality, that I(R) > 0. Let us denote 

(2.36) r = inf{r G (0, R) \ J(r) > 0}, 
and observe that, being I(r) increasing, we also have 

(2.37) I(r) = for every r < r' , 
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(2.38) I(r) > for every r,r' < r < R, 

(2.39) u(x) = M, for every x s.t. \x\ < r , 
for some Mel, 

Note that, by the maximum principle, if H(p) = for some p S (0, R), then u = in 
Bp and H(r) = for every r < p. Therefore H{r) > and N{r) is well-defined for every 
r G (r',fl]. 

Lemma 2.7 There exists a positive constant C , only depending on K and L, such that 
e p oN(r) is an increasing function on (r',R). 

Proof Let us recall that I(r) is continuously differentiable in (0, R) and that H(r) is 
absolutely continuous in (0,R), hence N is absolutely continuous in (r',R). If r £ (r',R), 
recalling Lemma 12.61 and using Schwarz inequality, we compute 

, , . N'(r) 1 I'(r) H'(r) 

(2.40) — — = — I — — > 

1 ; N(r) r I(r) H{r) ~ 

> J 9 B r L M Vu -^ 2 _ J 9 bA AVu - v ) _ c > _c 

where C > only depends on K and L. Therefore 

C 

(2.41) N'(r) + —N(r) > 0, for every r 6 (r', R), 

Po 



or, equivalently, 



d r^- 

(2.42) — (e >>oN(r)) > 0, for every r G (r',fl), 

which proves the lemma. □ 

Proof of Theorem 12.31 Let r%, r 2 , ?*3 be such that r' < n < r2 < r$ < R. By inequality 
(I2.25j) and by Lemma 12.71 we have 

x , H(r 2 ) f r2 H'(r) N , r 2 f r2 N(r) C , 

2.43 l g-p(=/ ~^rdr< n-1 log ^ + 2 / -^2 dr + - r 2 - n < 
ff(ri) i ri H(r) n J n r p 

< (n - 1) log - + e^ (r2 - ri) ^(r 2 ) log ^ + — (r 2 - n), 



(2 .44) log §^ . f "ft* > (n — 1) log — 4- 2 f ^* - ^ - > 
H(r 2 ) Jr 2 H{r) r 2 7 r2 r p 

> (n - 1) log — + e^ (r2 " r3) iV(r 2 ) l og Tl - —(r 3 - r 2 ), 
r2 r 2 po 
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where C > denotes a constant, only depending on K and L. Prom (|2.43[) and ()2.44j) . we 
have 



(2.45) log|M< (n _i)iog^ + 



where 



c_ 



(rs-n) / log H(r 2 ) , ,v . C(r 3 -r 2 )\ r 2 C 

H (n-l) + ; ft - 1o S — + — r 2 - n < 

^ cL (r3 _ ri) log^ iJ(r 3 ) C. , rf r ,- r ,) lo Sn c , 



log a ff(r 2 ) po "' log a Po 



< P iog^4+s, 

#(r 2 ) 



(2.46) p = e loe^^ ' 



(2.47) 5 = Cmax{ — ,1 Ul+p) 



Po 

and C > only depends on if and L. We have 

(2 48) H{r 2 ) ( H{r,)Y n 

(2 - 48) hW) ~ \hW)) e ' 



(2.49) (H(r 2 )) l+P < e B H( n )(H(r 3 )r, 

(2.50) H(r 2 ) < e Cmax {n' 1 }(H{r 1 ))^{H{r 3 ))^. 

Therefore, assuming the additional condition r\ > r 1 , inequality (|2.T|) follows with Q = 

exp (Cmax{|^,l}) and a = j^. 

It only remains to prove that r' = 0. Let us assume by contradiction that r' > 0. If we 
had M = in (|2.39[) . then H(r') = and, passing to the limit in (|2.7p as r\ — > r', it would 
follow H(r) = for every r 6 (r',R) and u = in .Br, leading to a contradiction with the 
definition (|2.36p of r' . Therefore in this case r' = 0. If M ^ in (|2.39p . let us consider the 
function v = u — M, which satisfies Vw = Vit in Br and v = in 5 r '. By applying the 
above arguments to v, we have again that r' = and the proof is complete. □ 



Remark 2.8 It is worth emphasizing once more the strength of the frequency function 
method just employed. In fact, besides the three-spheres inequality, it enables also to obtain 
a doubling inequality, and this was in fact the original purpose when the method was devised 
by Garofalo and Lin [32]. Let us outline here how the proof goes. For the sake of simplicity 
we stick to the normalized setting of Theorem 12.31 From the proof of Theorem 12.31 we have 
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obtained that r' = and, consequently, Lemma 12.71 applies for all r £ (0,R). Hence, by 
(f2T25j) and by Lemma O 



(2 . 51) log = t m ir < (n- i)b g 2+2 r a + £ < 

#(r) J r H{r) J r s po 



<( n -i)io g 2 +2 r e *" i -»^ (is+ ^< 

Jr 8 PO 



C R 

< (n - 1) log 2 + 2CN(R) log 2 + . 

Po 



That is 



(2.52) H(2r) < QH(r), for every < r < — , 

where Q > 1 only depends on K,L, ^ and on N(R). Consequently, by an integration and 
by the ellipticity bounds (|1.19f) . 

.2 ^ f „.2 - 



(2.53) / u 2 <Q u, for every < r < 

with a possibly different Q, but still depending on the same quantities. As is well-known, 
(|2.53p in turn implies the strong unique continuation property, because, by iteration on the 
radii r n = R2~ n , n = 1, 2, . . ., we readily arrive at 



f i 1 , r s '°bQ r o 
(2.54) / u z > — (— )io g a / u^, for every 

Js r <y R Jb r 



< r < -. 



3 The two-dimensional case 

In this section we shall deal with the two-dimensional case. We shall prove, by complex- 
analytic techniques, a three-spheres inequality for an elliptic equation in divergence form in 
pure principal part, see Theorem l3.111 The main difference with respect to higher dimensions 
is that in two dimensions we allow the coefficient matrix A to be discontinuous. Actually 
here we allow A to be also non-symmetric, however we stress that this hypothesis will be 
used in this section only. 

Throughout this section we shall identify z = x + \y £ C with points (x, y) £ M? and we 
shall fix A = A(z), z E M 2 , a real- valued 2x2 matrix such that its entries are measurable 
and it satisfies, for a given constant K > 1, the ellipticity condition 

(3.1) 1-^(1)1 £ > K~ '|£| 2 f ° r almost every z £ R2 ' for every £ G R2 - 

The constant will be referred to as the ellipticity constant of A. When A is symmetric 
this condition coincides with the one given in (11, 19ft . We also observe that (|3.1j) implies that 

|fljj(-2)| < A for every i, j = 1,2 and for almost every z € M 2 , 

A depending on i^T only. 
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In this section, f2 will as usual denote a bounded connected open set contained in M?. 
Letting / : 0, — > C be a complex- valued function, we shall make repeated use of the following 
notation for complex derivatives 



We denote by J = [° q 1 ] the counterclockwise rotation of \. 

Proposition 3.1 Under the above stated hypotheses, let u £ H} oc {ft) be a weak solution to 
the equation 

(3.2) div(^Vn) = 0, in Q, 

where Q, is a bounded simply connected open set in M 2 . 
Then there exists a function v £ iJ^ c (fi) that satisfies 

(3.3) Vf = JAVu, almost everywhere in O. 
Moreover, letting f = u + 'w, we have 

(3.4) fz = u f z + vf z , almost everywhere in S7, 

where fx and v are bounded, measurable, complex-valued coefficients satisfying 

(3.5) + \u\ < k < 1, almost everywhere in Q, 

A; &em(/ a constant depending on K only. 

On the other hand, if f = u + if £ Hl oc {£l,C) solves (|3.4p - (|3.5p . i/ien i/iere exists a 
real-valued 2x2 matrix A £ L°°(f2), satisfying (|3,ip urei/i a constant K depending upon k 
only, such that u is a weak solution to f|3.2|) . 



Proof By (j3. 3f) the vector JAVu is curl-free in the weak sense, thus v is well-defined as 
long as is simply connected, see [8, Theorem 2.1] for details. Then (|3.4p follows from (|3.3p 
with /x, given by 

_ «22 ~ «n ~ i(«i2 + Q21) 

Oll«22 — 012021 + «11 + a22 + 1 ' 

(3-6) 

_ Q12Q21 - «iia22 + 1 + i(ai2 ~ Q21) 
ana22 — 012021 + an + 022 + 1 

From these expressions and (|3.ip . one obtains, through elementary although lengthy com- 
putations, (|3.5p . In particular, is related to by the following formula 



K + VK 2 -1 - 1 
~ K + VK 2 -l + l' 

This estimate is sharp and it is proved in [9]. 
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Conversely, given the coefficients fx, v in (|3.4p satisfying (|3.5p one obtains (|3.2|) and (|3.3p 
with A given by 




(3.7) 



and the conclusion follows. It may be shown that the ellipticity constant of A is related to 
k by the following formula 

1 + k 



K < 



1 - k' 



Also this estimate is sharp, see again 



□ 



The function v appearing above is usually called the stream function associated to u. 
Notice that v is uniquely determined up to an additive constant and also that v is a weak 
solution to 



(3.8) div(BVv) = 0, in ft, 

where B = (detA)~ 1 A T , and where det(-) denotes the determinant and (-) T denotes the 
transpose. We notice that B still satisfies the ellipticity condition (|3.ip . with the same 
ellipticity constant. 

For any k, < k < 1, we say that a function / : ft — > C is a k-quasiconformal function 
in ft if it is an i?^ (ft, C)-solution to the Beltrami-type equation (|3.4p - (|3.5p . A univalent 
solution to (I3.4p - (l3.5p is said a k-quasiconformal mapping. A function / is a quasiconformal 
function, respectively mapping, if it is a k-quasiconformal function, respectively mapping, 
for some k, < k < 1. Concerning quasiconformal functions, their properties and charac- 
terizations we refer to [66]. 

Finally we observe that when A is the identity matrix, then u is harmonic and v is 
simply its harmonic conjugate, thus / = u + iv is a holomorphic, or conformal, function, 
corresponding to k = 0. Such a special case deserves our attention, especially when ft = B\. 
We have the following result which links properties of a harmonic function u to properties 
of its harmonic conjugate. 

Proposition 3.2 Let u be a harmonic function in B\, and let v be its harmonic conjugate 
normalized in such a way that v(0) = 0. 

Let us assume that \u(z)\ < E for every \z\ < 1. Then we have 

, s | / \ I 2E , (\ + \z\\ 

(3.9) |«(z)|< —log^UJ. 

If u is Holder continuous on B\, with Holder exponent /3, < (3 < 1, then v is also 
Holder continuous on B\, with Holder exponent j3 and 

W v Wc0(B[) - C o\\ u \\ct>(B[y 

Co depending on f3 only. 
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Proof If u is continuous up to the boundary of B\, then 

f(z) = u(z) + iv(z) = tt-t / T^— U (€)-T, for ever y M < !• 

Such a representation is often referred to as Schwarz's formula or as Poisson represen- 
tation. The inequality (j3.9f) follows immediately and it is sharp, take for instance u = 
9 ( - log I y±| 11 on 5i. The second result is due to Priwaloff , see [87] . □ 



Properties of harmonic functions and their conjugates, or of holomorphic functions, 
may be transferred to solutions to elliptic equations in divergence form and their stream 
functions, or to quasiconformal functions, through the following important representation 
theorem, a proof of which is due to Bers and Nirenberg [18], see also [17]. 

Theorem 3.3 Let £1 be a connected open set contained in B\ and let f £ H} (Q, C) solve 
(|3.4p where \i, v satisfy (|3.5p . 

Then there exist a k- quasiconformal mapping x from B\ onto itself and a holomorphic 
function F on x(^) such that 

(3.10) f = Fo X . 

We may choose x such that x(0) = 0- Moreover we have that the function \ an d its 
inverse x~ 1 satisfy the following Holder continuity properties 

(3.11) |x(x) -x(y)\ < C\\x-yf, for any x,y G B 1 
and 

(3.12) \x' 1 {x)-X~ 1 {y)\<Ci\x-y\ 13 , for any x,y G B u 
where C\ and (3, < (3 < 1, depend upon k only. 

Proof See [18, page 116] and [17]. □ 

Through the use of quasiconformal functions it is also possible to extend the classical 
method of harmonic measure in order to obtain a Holder stability estimate in the interior for 
Cauchy problems for solutions to Beltrami- type equations (|3.4|) - (j3.5p . In turn, this provides 
a Holder stability estimate in the interior for Cauchy problems for solutions to (j3.2[) . 

We denote by Ca the differential operator 

(3.13) C A u = div(AVu). 

We now introduce the notion of /^-harmonic measure w of a subset S of d£l. Roughly 
speaking uj is the solution to the Dirichlet problem 

(div(ylVLj) = 0, in ft, 
w = l, on S, 

oj = 0, ondn\ S. 

However, being the set S arbitrary, and since the Dirichlet data in (|3.14p may be discontin- 
uous, a more careful definition is needed. Let us recall some notions from potential theory, 
see for instance [41]. 



26 



Definition 3.4 A function u : 0, R u {+00} is called CA-superharmonic in Q if 

(i) u is lower semicontinuous; 

(ii) u ^ +00 in any connected component of £1; 

(iii) for any open set fii CC and any h £ C(f2i), such that £,4/1 = in the weak sense 
in Oi, if u > h on dSli then u > /t in 

A function u is CA-subharmonic in $7 if — u is £^-superharmonic in f2. 

Definition 3.5 Let E be a subset of d£l and let xs t> e its characteristic function. We 
define Us as the class of the /^-super harmonic functions u in fl such that u > and 
liminf^^y > Xs(y) f° r an Y y £ <9f2. 

We define the CA-harmonic measure of E with respect to O as the upper Perron solution 
with respect to X£> that is 

lu(z) = w(E, fi, >Ca; 2) = inf{n(z) | u £ ^s}, for any 2 6 O. 

We observe that we have < uo < 1 everywhere. In order to have that < to < 1 in Q, we 
need to guarantee that to is not identically equal to or to 1. A simple sufficient condition 
is the following, [41, Theorem 11.6]. We recall that a continuum is a closed connected set 
with at least two points. There exist two continua, 71 and 72, such that 71 is contained in 
E and has positive distance from dVt \ E and 72 is contained in d£l \ E and has positive 
distance from E. For a more thorough analysis of this question we refer to [41]. 

Lemma 3.6 Let f £ H} (O, C) be a non-identically zero solution to (|3.4p - (|3.5p . Then there 
exists a real-valued 2x2 matrix A\ £ L°°(£l) such that A\ satisfies (|3.ip with a constant 
K depending on k only, and cp = log |/| is CA 1 subharmonic. 

Proof We denote 

_ f [i, _ where f z = 0, 
/iX 1 M+ vfz/fz, where / 2 / 0. 

Then we have that \ni\ < k almost everywhere in Q and / satisfies 

h = Ml/z> almost everywhere in Q. 

Let z be a point in Q such that /(z) 7^ 0. Locally, on a neighbourhood of z, we can 
define the function = log / where log is any possible determination of the logarithm in 
the complex plane. In this neighbourhood <!> also satisfies 

(3.15) *z = A*l**- 

Then by Proposition 13. II the function tp = log |/| = 3? log / locally satisfies 

(3.16) div(AiV<£) = 

where the matrix A\ is given by (|3.7|) with \i and v replaced by \i\ and 0, respectively. 
Note incidentally that A\ is symmetric, that is A\ = Aj, det^i = 1, and hence B\ = 
(det^i)" 1 ^ = A x . 
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We remark that we can define ip = log|/| globally as an Hh c (^i) function, where 
Oi = {z G f2 | /(z) 7^ 0}, hence using a partition of unity it is easy to show that (|3.16p 
holds weakly in Oi. 

Clearly, by Theorem 13.31 the se t {z £ £1 \ f(z) = 0} consists of isolated points and f 
goes uniformly to — oo as z converges to an element of such a set. 

Using this remark and the maximum principle we can prove in an elementary way that 
if = log | /| is £Ai-subharmonic. □ 

We observe that if k = 0, that is / is holomorphic, then A\ in the previous lemma is 
the identity matrix. By the use of the Ca 1 -harmonic measure we obtain a Holder stability 
estimate in the interior for Cauchy problems for a Beltrami- type equation (I3,4p - (l3,5p , 

Theorem 3.7 Let E be a subset ofdSl. Let f G #/ oc (ft,C) solve (|3T3]) . 
V H/llL°°(n) — E and we have that, given r\ > ; 

(3.17) limsup |/(x)| < rj, for any y G E, 

x^y 

then for any z G f2 the following estimate holds 

(3.18) |/(z)| <-E 1 ~ w( *V ( *\ 

where uj = w(E, Q,/^) is i/ie £,A 1 -harmonic measure o/E mt/i respect to O and £/ie matrix 
A\ is defined as in Lemma \3.b\ 

Proof This result was obtained in [10, Theorem 4.5]. We sketch a proof, for the sake of 
completeness. It is evident that < r] < E. If n = E the proof is trivial. We consider 
then the case < rj < E. By Lemma 13.61 the function f = log |/| is £/ii-subharmonic. Let 
to = cj(E, £1, Ca x ) be the /^-harmonic measure of E with respect to O. 
Let us denote 

log (Jf 

<P 

2L 



It is easy to see that f belongs to the upper class U-£. Hence for any z G 0, we have 
w(z) < and consequently 

(3.19) p(z) < logfaM*) + log(£)(l - lo{z)). 

Thus the conclusion follows. □ 



We observe that, in view of Proposition 13.11 the above Theorem 13.71 could be restated 
in terms of a Cauchy problem for an elliptic equation like (|3.2p . 

Theorem 3.8 Let E be a subset of dd. Let u G be a solution to (|3.2p . Let us 

assume that there exists a single-valued stream function v in Q and let f = u + iv. We 
a-priori assume that < E. 

We assume that E C dVL is a connected open Lipschitz arc of dQ and as usual we denote 
g = u|s and ip = AS/u ■ the Cauchy data of u on E. We suppose that g G C°(E) and 
ip G L X (E). 
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//, given n > 0, 

NU«(E) + MUi(E) < V, 
then for any z£l! the following estimate holds 

(3.20) \u(z)\ < E^WrfW, 

where u = w(E, Q,Cai) is the C ^-harmonic measure o/E with respect to 0, and the matrix 
A\ is defined as in Lemma \3.b] with fj, and v given by (|3.6|) . 

Proof By our assumptions on g and ip, and classical regularity estimates, we have that / 
is continuous at any point y G E. 

We normalize v in such a way that v = on some point of E, and we infer that 
||/IU°°(£) — V and, in view of the regularity of /, we obtain that (|3.17p holds. Therefore the 
conclusion is an immediate consequence of Theorem 13.71 □ 

It is evident that in order to make such estimates practically useful, it is necessary to 
provide a positive lower bound on uj and this is usually obtained by using the Harnack 
inequality for positive solutions to elliptic equations in divergence form, [34, Theorem 8.20], 
see for details [10]. 

Our principal aim here is however to concentrate on three-spheres inequalities, first for 
holomorphic and quasiconformal functions, then for harmonic functions and solutions to 
elliptic equations in divergence form. 

Theorem 3.9 Let R > and < k < 1, and let f be a k- quasiconformal function on Br. 
Then for every r\, r<i, r^, with < r\ < r<i < r$ < R, there exists a constant a, < a < 1, 
depending on k, ^ and ^ only, such that 

(3-2i) ll/lU-(B r2 )<ii/ii2oc (Bri )ll/lliq Br3 )- 

Proof We may use the previous Theorem [321 by setting 0, = B r3 \B ri and E = dB ri C d£l. 
Then 

a = inf u(H,CI,Cai)(z)- 

zeB r2 \B ri 

For example, when k = 0, by this technique we may recover the three-circles theorem 
for holomorphic functions by Hadamard. In fact, in this case A\ is the identity matrix and 
by explicit computations in radial coordinates we obtain 

log(^) 

a = i — when k = 0. 
log(^) 

In order to obtain an essentially optimal value of a in the general case, that is when 
< k < 1, we use the representation theorem, Theorem 13.31 A related argument was used 
already in [7, Proposition 1] . Up to a linear change of variables, we may temporarily assume 
that r3 = 1. We apply Theorem 13.31 then / = F o x, where F is conformal on B\ and 
X is a quasiconformal mapping between B\ and itself. Let C\ and j3, < /3 < 1, be the 
constants appearing in (13.111) and ()3.12j) . which characterize the Holder continuity of the 
quasiconformal mapping % and of its inverse. Then there exist < r i < ?2 < 1 such that 
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-Bf-L C x(B ri ) and x(B r2 ) C B^ 2 or, as is the same, x{B\ \ Br-i) D B\ \ B^ 2 . We apply the 
three-circles theorem to F, and we obtain that (|3.2ip holds with 



Q = 

where we may take 



(3.22) a = when < k < 1, 

l°g(^) 



7"3 



(3 ' 23) ?1 = feJ md ' 2 = 1 \(h 

C\ and /3, < j3 < 1, depending on k only. The above choices of r% and ?2 serve our purposes 
for any r% > 0. □ 

We now turn our attention to the second order elliptic equation. We shall prove a three- 
spheres inequality first for L°°-norms and then for L -norms. Let R > and let u G H} oc {B R ) 
solve 

(3.24) div(,4Vu) = 0, in B R . 

Theorem 3.10 Under the above stated hypotheses, for every r%, r^, r%, with < r\ < T2 < 
i~3 < R, there exist constants Q > 1 and a, < a < 1, depending on K , ^ and £| on/?/, 

(3-25) IMU~(B r2 ) <Qll«ll£-(B Pl )||tt|li« a (Bfa) - 

Proof We temporarily assume that r% = 1 and also ^4 = Id, that is ?i harmonic in B\. We 
let ?; be its harmonic conjugate normalized in such a way that ?;(0) = 0. By (|3.9p . we may 
bound the values of v in terms of those on u. For example, we have 

ll u IU°°(Bg_) < C2||w|U°°(B ri ) and ||?^||x°°(g r2+ i) < Cs (Si) 



with 



2 . „ 2 /3+ 



12. 



(3.26) C 2 = -log3 and C 3 = - log ,.. 

7T 7T \ 1 / 

Therefore, we may use the Hadamard three-circles theorem and conclude that for u harmonic 
in B R we have 

(3.27) [|u||l»(b„) < ((C 3 + l)||«|Uoc (Br3) ) 1 -«((C 2 + l)|MU«, (Bpi) ) Q 
with a given by 



(3.28) a 



log (s + £) 



For the general case of solutions to elliptic equations in divergence form, we use the same 
argument as in Theorem 13.91 By Proposition 13.11 and Theorem 13.31 we have that u = U o\ 
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where U is a harmonic function in B\ and x ls a quasiconformal mapping between B\ and 
itself. Then we apply (|3.27p to U and we conclude that (13.251) is satisfied in the following 
form. It holds (|3.27[) under conditions (|3.26p . (|3.28[) . with ^ and — replaced by T\ and ¥2, 
respectively, ¥\ and ¥2 given by (|3.23p . □ 

We remark that for harmonic functions, instead of (I3.27h . we might have used a result 
by Korevaar and Meyers, [56]. They proved that for u harmonic in Br a three-spheres 
inequality holds with Q = 1, that is 



(3.29) IM|ioo (Bf , 2 ) < \\u 



il-5 



r 2 , - i\L°°(B r3 )\\ U \\L°°(B ri )> 

where 

(3.30) a = a ( — , — 

We observe that their result is actually valid in any dimension n > 2. Let us also recall an 
interesting related result by Petrosyan [85]. 

As a corollary to Theorem 13. 1 II we obtain the corresponding three-spheres inequality in 
L 2 -norms as well. 

Theorem 3.11 (Three-spheres inequality — two-dimensional case) Under the pre- 
viously stated hypotheses, for every r\, r2, r%, with < r\ < r2 < r% < R, there exist 
constants Q > 1 and a, < a < 1, depending on K , — and ^ only, such that 

( 3 - 31 ) \\u\\mB r2 ) < QhW^^JluWl^^y 

Proof We use the local boundedness estimate for weak solutions to (|3.24p . see for instance 
[34, Theorem 8.17], telling that for a constant C depending on K only, and for every s, r, 
< s < r < R, we have 

II II < d P° II II 

\\u\\l°°(b s ) < — -IMIl^b,)- 

We apply this bound with the following choices of s, r 

r2 + r 3 

s = . 

2 

and also 



^3, 



s = —, r = n, 

in combination with Theorem 13.101 with radii < ^ < r2 < r2 + r3 . Hence, by the trivial 
estimate _ 

II II ^ V^ r 2n I, 

\\ u \\L2{B r2 ) < —j^-\m\L°°(Br 2 )> 

we conclude that 

\H\l\b T2 ) < V^QC (j^^MlHb,,^ (^-\\u\\ L 2 (Bri ^j 
where Q and a are the constants appearing in (I3T25|) related to K, 2a an d r -^. □ 
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4 The three-spheres inequality for the complete equation 

Here and in the rest of the paper we consider the elliptic operator div(^4V-) +c • where A and 
c satisfy the conditions (|l,19p . (|1.20p . and (jl.2ip . respectively, stated in the Introduction. 
We recall that the Lipschitz condition (|1.20p is invoked only when n > 3. 
We start by considering solutions u to the homogeneous equation 

(4.1) div (AVu) + cu = 0, in B R . 

Theorem 4.1 (Three-spheres inequality — equation with zero order term) If the 

above stated hypotheses hold, there exists Cq, < Cq < 1, only depending on K , L and n, 
such that, setting 

(4.2) Ro = mm{R,C QPo }, 
for every t\, r2, r^, with < r\ < ri < jf^ < r^ < Rq, 



(4-3) \\u\\m Br2 ) < Q\W\ 



L*(B ri )\\ u \\L2(B r3 y 



il-a 



where Q > 1 and a, < a < 1 only depend on K, L, k, max ll, £a, 
The thesis is an immediate consequence of the following two lemmas. 

Lemma 4.2 For every 5 > there exists Cq, < Cq < 1, only depending on K , L, k and 
5 such that, denoting Rq = min{R, CqPq}, there exists a positive solution w £ C 1 (Bji ) to 

(4.4) div (AVw) + cw = 0, in B Ro , 
such that 

(4.5) y^<w<\ + 5\ in Br,, 
and when n > 3 

(4.6) \Vw\<— , in Br,. 

Po 

Proof There exists C\ > only depending on K, k, such that the operator div (^4V •) + c • 
is coercive on Hq(B t ) for every r < R\, with R\ = min{i?, C\Pq}, see for instance [34, 
Lemma 8.4]. Let us pick the unique solution w to 

div (AVw) + cw = 0, in B r , 
w = 1, on dB r . 



div (AVz) = -c(l + z), in B r , 
z = 0, on dB r . 
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By standard a-priori bounds in L°°, see [34, Theorem 8.16], we have 

(4-9) IWIi°°(B r ) < Ck— (1 + \\z\\ L *o {Br) ), 

Po 

where C only depends on K. Hence there exists C% < C\ only depending on K, k, such 
that for every r < R2, with R2 = min{i?, C2P0} we have 

(4.10) INIU-(B r ) < C K~2' 

Po 

Next, when n > 3, applying to (|4.8|) a global Schauder type estimate, see [34, Theorem 8.33], 
one obtains, for every r < R2, 

r 

( 4 -ll) ||Vz[|ioo( Br ) < Ck— (1 + ||z||L°o( Sr )), 

where C only depends on if, L. Thus, we may find Co < C2 only depending on K, k, L 
and 5, such that, denoting i?o = min{i?, Copo}, we have 

( 4 - 12 ) \\4l^(B Ro ) < 

and when n > 3 

(4-13) l|V^||^(^ ) < 



Therefore (|4.5p and (|4.6p follow immediately. 

Observe that, for the sake of simplicity and with no loss of generality, we can also assume 
Co < 1. □ 



Lemma 4.3 Let Rq be the quantity introduced in Lemma \4--S\ when 5 is chosen to be 5 = 1 
when n = 2 and 5 = min{-^, 1} when n > 3. Then u can be factored in Br as 

u = wv, 

where w is the function constructed in Lemma \4-%\ and v solves 

(4.14) div(lVv)=0, inB Ro , 
where the matrix A is given by 

(4.15) A = w 2 A, 

and it satisfies 

1 ~ 

(4.16) ttjI^I 2 < A(x)£ ■ £ < 4if|£| 2 , for almost every x 6 -Br , for every £ £ M n . 
and when n > 3 

~ ~ 8-L 

(4.17) \A(x)-A(y)\< — \x-y\, for every x,y £ B Ro . 

Po 
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Proof The proof is straightforward. □ 

Let us now consider the inhomogeneous elliptic equation 
(4.18) div (AVu) + cu = f + div F, in B R , 

where / G L 2 (R n ) and F G L 2 (K n ;R n ) satisfy (f!~29l) . 

Theorem 4.4 (Three-spheres inequality — complete equation) If the the above sta- 
ted hypotheses hold, for every r\, r<i, r%, with < r\ < ri < 4^ < r% < Rq, 



l-a 



( 4 -!9) \\ u \\LZ(B r2 ) < Q (ll«llL2(B ri ) + s) (|MU 2 (B r3 ) + £ 

where Q > 1 and a, < a < 1, on/y depend on K, L, n, max|^,l|, and i?o is 

given by (|4.2p . 



Proof Let us consider the unique solution uo to 

div (^4V«o) + cuo = / + div F, in £?# , 

uq = 0, on <9-Br . 



(4.20) 



We have that 

(4-21) II u o||l 2 (B Bo ) < C(-RoII/IIl 2 (M' 1 ) + ^o|I^IU 2 (R";K«)) 

with C only depending on K. Noticing that u — uq satisfies the hypotheses of Theorem 14. 1\ 
the thesis follows immediately. □ 

5 Propagation of smallness 

For every G C W 1 and for every h > we shall denote 

(5.1) G h = {x G G | dist(x, 0G) > h}, 

(5.2) G* 1 = {x G R" dist(x,G) < h}. 

Theorem 5.1 (Propagation of smallness in the interior) Let Q be a bounded con- 
nected open set in R n , and let B To (xq) C be a fixed ball. Let Cq be as in the thesis 
of Theorem \4-l\ Let h, < h < min{2CopO) oe fixed and let G C be a connected open 
set such that dist(G, dU) > h and -Bro(xo) C G. 
Let u G Hi oc (Q) be a solution to the equation 

(5.3) div (AVu) + cu = f + div F, in tt, 
where f and F satisfy (|1.29p . Let us assume that 

(5-4) h\\L2(B ro (x )) < V, 
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(5-5) IMlL 2 (n) < 

for given ij > 0, Eq > 0. We have 
(5.6) [Ml 1,2 (G) < C (77 + r)' 1 (Ed + -) 



(5.7) C = Ci 

and 



1 



7i"~ 



c 2 |ni 

(5.8) <5 > aT 

u>it/i Ci > and a, < a < 1, only depending on K , L and k and C2 only depending on 
K. 

Proof We shall need uniform three-spheres inequalities in a domain slightly larger than G. 

For instance, for any x G G2, we have Bh(x) C O. Therefore we can apply the three-spheres 

2 

inequality (|4.19p to spheres centered at x choosing R = h/2. Hence Rq = min{^, C0P0} and, 
by our choice in the assumptions, i?o = f • Moreover, recalling that Co < 1, we also have 

max {S L ' 1 } = L 

Next we can fix radii n, r$ as follows 

/i r 3 /i 1 /» 

( 5 - 9 ) r 3 = -, r 2 = ^7 = 77-7, ri = -r 2 



2' 5if 10K' 3 30K 

With such a choice the inequality (|4.19p applies with Q > 1 and a, < a < 1, only 
depending on K, L and k. 

Let us consider the set G ri as defined in (|5.2p . We have that G ri is a connected open 
set containing G such that dist(G ri , <9Q) > /i — rx > r 3 = |. For every y 6 G"" 1 , there exists 
a continuous path 7 : [0.1] — > G"" 1 such that 7(0) = xq, 7(1) = y. Let us define = to < 
t% < . . . < tN = 1, according to the following rule. We set t^+i = max{i | [7(f) — x k \ = 2ri} 
if \xk — y\ > 2n, otherwise we stop the process and set N = k + 1, izv = 1. Let = 7(ifc). 
The balls B ri (xk) are pairwise disjoint for fc = 0, ... , JV — 1 and |x&-fi — x^l = 2ri. Since 
i"2 = 3ri we have that B ri (xk+i) C B r2 (xk) and therefore, by the three-spheres inequality 

(iron . 

(5-10) IMU 2 (B ri (*fc+i)) +£ - ^ (IMIl^ +e) (^o+e) 1 "", 

for k = 0, . . . , N — 1, where Q > 1 and a, < a < 1, only depend on K, L and k. 
Denoting 

( Xfc )) +£ 

mfc = E^Ts ' 

we then have 

(5.11) m k+1 < Qm%, for k = 0, . . . ,N - 1, 
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(5.12) m N < Qm s , 

where Q = Q 1 + a +—+ a 1 an d § = a N . Hence we have obtained 

(5-13) h\\L*(B ri (y)) < Q (\\u\\L*(B ri (x )) +e) (-Eo + e) 1_<5 - 

Now, 1 + a + . . . + a^ -1 < Since B Tl (xq),. . . ,B ri (xn-i) are pairwise disjoint, we have 

that N < < g ^lp^ , with C2 only depending on i-T. Hence, recalling that Q > 1, we 

compute 

(5.14) g<Q^, 



(5.15) 5 > a >>™ . 

Let us tessellate M. n with internally non-overlapping closed cubes of side I = ^= and let 
Qj, j = 1, . . . , J, be those cubes which intersect G. Clearly, any such cube is contained in 

n 

a ball of radius r\ and center Wj G G V1 and J < n 2 n r n ■ Therefore, from (I5.13|) . we have 

(5.16) / < W «* < £ / «* ^ J Q 2 Po (MlHb^o)) +ey 5 (E + e)^- s l 

JG ~{JQj ~[JB ri {w 3 ) v ' 

Hence the thesis immediately follows. □ 



Remark 5.2 It is important at this stage to emphasize that the above Theorem 15.11 on 
propagation of smallness in the interior enables us to generalize the three-spheres inequal- 
ity of Theorem 14.41 by removing the limitations on the radii that were present there. It 
should also be mentioned however that by this approach it does not seem possible to obtain 
optimal estimations of the constants C, a appearing in the inequality (|1.42|) . This optimiza- 
tion problem for such a general version of the three-spheres inequality still mantains some 
unanswered aspects. 

Proof of Theorem II. 101 It follows by applying Theorem l5-H with B ro (xo) = B ri , G = B r2 , 



The dependence of C, a on the quantities stated in Theorem 11.101 is straightforward, 
although somewhat lengthy. □ 



Theorem 5.3 (Global propagation of smallness) Let f2 be a bounded connected open 
set in W 1 with boundary d£l of Lipschitz class with constants po, Mq. Let u £ H 1 ^) be a 
solution to the equation 

(5.17) div (AVu) + cu = f + div F, in SI, 

where f and F satisfy (jl,29p . Let B To {xq) C SI and let us assume that 
(5-18) \\u\\L*{B ro [ xo )) < V, 
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(5-19) IMIffipi) < E, 

for given n > 0, E 1 > 0. We have 

(5-20) || n || i2(n) <(£ + e )^|±£^ 

u;aere 

C_ 



(5-21) ^)<- — rni, fort<l. 



where C > and < p, < 1, only depend on K, L, k, Mq, ^ and jpr. 

Remark 5.4 As it will be evident from the proof, rather than the bound (|5.19p we shall 
actually use a weaker one. In fact we might replace (|5.19p by the assumption that there 
exists a p > 2 such that 

(5-22) Nl£P(n) < E. 

In this case the constants C, p would also depend on such an exponent p > 2. We have 
chosen to formulate Theorem 15.31 (and also Theorem 1 1 .9j) in terms of an a-priori i? 1 -bound 
because in many applied settings such a bound has a clearer physical interpretation. 

We premise the following proposition. 

Proposition 5.5 Let Q be a bounded connected open set in M. n with boundary d£l of Lips- 
chitz class with constants po, Mq. There exists ho > 0, only depending on pq, Mq, such that 
O/j is connected for every h < h$. 

In order to prove the above proposition, let us introduce the following lemma. We recall 
our assumptions. Q is a bounded connected open set in M n with boundary dO, of Lipschitz 
class with constants po, Mq. Let P £ d£l. According to Definition ll.il up to a rigid motion, 
we have P = and 

finr_Po_ on (P) = {x = (x',x n ) eVpo_ i x n > z(x ; )}, 

where Z is a Lipschitz function on B' PQ C M n_1 satisfying 

Z(0) = 0, 
\\ Z \\ ( \< MqPq. 

We denote 

(5.23) d{x) = dist(x,dfi), x e n, 

(5.24) d(x) = x„ - Z(x'), xefinr 
and also 

(5.25) do- P ° 



2M 
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Lemma 5.6 Under the above stated assumptions, if we assume 

(5.26) x £ T_n_ nPin, d(x) < d , 

2M '2 

then we have 

(5.27) d{x) < d{x) < (l + 1 + M^j d{x). 

Proof Let x £ T_p^_ po n Q be such that d(x) < do. Let z £ d£l be such that \z — x\ = d(x). 

2M '2 

Recalling that Mq > 1, we have that 

I 'I s I ' 'I i I 'I / At \ i P ° S P0 

\z\<\z -x| + |x|<d(x) + ^<— , 

|Zn| < Fn - + \x n \ < d{x) + — < — — + — < p . 

Z ZlViQ I 

Therefore z £ Tpo_ po n d£l, so that z = (z' , Z(z')). 
Let y = (x' , Z(x')). We have that 



d[x) = \x — z\ < \x — y\ = d(x), 

d(x) = \x — y\ < \x — z\ + \z — y\ = \x — z\ + (\x — z'\ 2 + \Z{x) — Z(z')\ 2 )^ < 

<\x-z\ + yJl + M 2 \x - z'\ < M + yj 1 + M^j d{x). 

□ 

Proof of Proposition 15.51 Let ho = — ^ ^ =^ . We trivially have that /i < g^j. 

Let h < h$. Let us first prove that $7^ is connected. Given x, y £ O^, let 7 be a path in fi 
connecting x and y. If 7 is not contained in Q^, let us modify 7 to a new path 7 connecting 
x and y and contained in fi^. Let C = {^(t) \ d(j(t)) < h}. For every z £ C, there 
exists £ 2 S 90 such that 1.2 — £J < /1 so that 2 £ T _po_ po(£ z ), where T_po_ po (£ 2 ) is defined 

8A/ ' 8 8A/ ' 8 

according to the Definition 1 1. 11 Therefore the cylinders T _po_ po(£, z ) provide an open covering 

8M ' 8 

of C and, since C is compact, there exist zi,...,zn such that C C U^-,r_po_ m(£ z ). Let 

8M ' 8 

us denote for simplicity £j = £ 2i , T* = T_po_ po (£i). Note that each such cylinder is possibly 

8Mq ' 8 

oriented with respect to a different coordinate system. For every point ^(t) £ C n T 1 let us 
replace j(t) with 7(f) = (V(t), Z(i(t)) + (1 + y'l + M 2 )/i). We can apply LemmaESlto 
7(f), obtaining that 7„(i) - Z(i(t)) = d(j(t)) < (1 + ^1 + M$)h. Recalling that M > 1, 
we have that 

\l'(t) - (fc)'| = |Y(t) - &)'| < 3T7- < 



8M 2M 



|7n(t) " £| < l7n(t) " 7«(*)l + l7»(t) - Cnl < (l + ^1 



1 + Mo 2 )/ i + ^<^. 



2 
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Therefore j(t) E Tj^ n (^) Dtt and d(j(t)) < py(i) - (Y(t), Z(i (t))\ < (1 + ^1 + M$)h < 

2AIq ' 2 

do. We can therefore apply Lemma [5U1 to j(t) and, by (|5.27j) . 



i + yr+ 

The connected components of 7 _1 (C n T 1 ) are closed intervals I a . In order to glue together 
7 and 7, let us add, for each endpoint t of any I a and for any t E 7~ 1 (C n OF 1 ), the closed 
segment joining {-/'(t), 7„(i)) and (7'^), Z(Y(t)) + (1 + i/l + M$)h). Now we repeat the 
above arguments to the path so modified and to the cylinder T 2 , and so on for a finite 
number of steps. 

We can now conclude proving that Qh is connected for every h < h^. Let x, y E £lh 
and let /i' such that h < h! < min{d(x), <i(y), /io}- We have just shown that O^/ is path 
connected, hence there exists a path in fi^/ C £lh joining x and y. □ 



Lemma 5.7 Let Q, be a bounded connected open set in W 1 with boundary dVt of Lipschitz 
class with constants po, Mq. There exists a constant C > 0, only depending on Mq, such 
that 

(5.28) n h ) < c\n\— . 

Po 

Proof We refer to [11, (A. 3)] for a detailed proof. In fact this result was proved there under 
stronger regularity assumptions, but the same arguments apply also under our Lipschitz 
regularity assumption. □ 

Proof of Theorem [53] Let hi = min / Z&. 2CqOq\, where h = ^ «r has 

been introduced in Proposition 15.51 We have that is connected, it contains B^o (xq) and 
also we can apply Theorem 15. 1| with G = f2/ ll . That is 

(5.29) \\u\\ L 2 {nhi) <c(\\u\\ L 2 {B ^ {xo)) +e} (\\u\\ L2{n) + e) 1 " 5 < C(E + e) (f^) » 

where C > and 8, < S < 1, only depend on iT, L, k, Mq, ^, and 

Let r E (0, /ii). Given any point x G Cl r \£lh 1 , that is such that r < c£(x) < /ii, let £ G 90 

be such that d = d{x) = |£ — x\. 

According to the Definition ll.l| up to a rigid motion, fl n Tpo_ po (£) = {y = (y',Un) G 

r pq I y n > Z(y')}, with Z a Lipschitz function on B'p^ C M n_1 satisfying Z(0) = 0, 

||Z|| c o,i < M qPq. 

By (pT2?l) . we have that r < d < d(x) = x n - Z(x') < (1 + ^/l + M^)d. Let 10 = 
(x',Z(x')) G 90. We have that \x'\ < d, \Z(x')\ < M d. 

Let us translate the origin of the coordinate system from £ to w. 

Let *° = i^/SLi ^ - M ° hi ^ s ° = 3 ( i+v^+X 2 " % ) and yo = (0 ' ^ o) ' We have that 

4s + hi= to , 
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so that 



to + 4s + hi < po ~ M Q d, 



B Aso ( yo ) c n hl nr_P±_ diPO _ Mod (w) c n hl nr^ )PO (0- 



Let C M = {y = (y',y n ) I 2/n > M|?/|} be the open cone with vertex at u; and tangent to 
the ball B so (yo), that is M = with sini? = ^. Let us denote s = sin$, q = Let 

si = gso, ti = gto and yi = (0,ii). Then the ball B si (yi) is tangent to the cone C w and to 
the ball B So (y ). By induction, setting s k = gs fc _i = q k s , t k = qt k _i = q k t , y k = (0, t k ), for 
k > 2, we have that the ball B Sk (y k ) is tangent to the cone C w and to the ball B Skl (y k _i). 
Let 

\\ u \\L2(B 3k (y k ))+£ 

mk = eTS • 

By applying (jl,42p to the balls of center y k and radii s k , 3s k , 4s k , and noticing that 
B Sk+1 (yk+i) C B 3sk (y k ), we have 

(5.30) Ttfc+i < C m fc > f° r every A; = 0, 1, ... , 

(5.31) m N < C l+a+ - +aN ~ 1 mf , for every N = 1, 2, . . . , 
where C > and a, < a < 1, only depend on K, L and k. By (15. 29ft we have 

(5.32) m <C '' ; 



E + e 

so that 



(5.33) m^v < C 



rj + e 



where a G (0, 1) only depends on K, L and k, whereas C > and (5, < 5 < 1, only depend 
on K, L, k, M , ^ and j|. 

In the new coordinate system centered at w, x = (0,d(x)) and, by (|5.27p and by the 
choice of to, r < d(x) < (1 + ^1 + M$)d < (1 + \A + Mq)^i < to- Since x O ftl , we have 
that rr B SQ (yo) U B si (yi) C i?3 SO (yo) C Hence there exists JV 6 N, JV > 2, such that 
x G B SN (y]y), x B SN1 (yjv-i)- It follows that d(a?) < tAr_i so that g^ -1 > ^. 

By (021) and by (15331) . we have 

II II / \ 'vryN- 1 

(5.34) NU^^)) < < / ^ 

where C > and 7 = a5 only depend on ET, L, k, Mq, ^ and Since a" -1 = q^ N > 

log a 
log 13 



^ 1 , and noticing that < 1, we have 

(5-35) ihiw^^)) < c ( E + £ ) (fri) nt " 
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with D = g and C > and 7 G (0, 1) only depending on K, L, k, Mq, £& and (§. 
Moreover sn = Q so > ^t— . Since B2 SN (x) C B3 Sn (i/n), we have that 



log 9 

to 

r \0 



(5-36) \\u\\ LHB2m{x)) <C(E + s)^^^j 



with C > and 7 G (0, 1) only depending on K, L, k, Mq, ^ and j^f. Let us tessellate W 1 
with closed cubes of side I = and let Qj, j = 1, . . . , J, be those cubes which intersect 
Q r \ f^. Clearly, any such cube is contained in a ball of radius and center such 
that |wj — Xj| < for some xj G Sl r \ $7^. By the above arguments, there exists JV£N, 
N > 2, such that |5.36p holds for x = Xj, with s^r > Therefore B SN (ujj) C B2 SN [xj) 

and, by (|5.36p and by the trivial estimate J < ^ 2 qrso)" ' we nave that 
(5.37) 

/ " 2 sE/ » 2 sE/ - 2 £ <W<* + -> 2 (f ) ~" 27fe) ° ■ 



where C > and 7 G (0, 1) only depend on K, L, k, Mq, ^ and ^. By (IQ9|) and by 
(|5.37p . we arrive at 



(5-38) \\ u \\ L2{nr) < C (E + e)[- 



r \ 2 / rj + e \ 7 V*o 



t / \E + e 



with C > and 7 G (0,1) only depending on K, L, k, Mq, ^ and j^f. Let us choose 

p = r^rg, for n > 2, whereas for n = 2, let us choose as p any number satisfying p > 2. From 
Holder inequality, Sobolev embedding theorem and by Lemma [ST71 we deduce 

. , , /in\n r |\5-* 

(5.39) \\u\\ L 2 in \ nr) < ^ — — — 1 ||«||jy>(n) < 

with C only depending on Mo and Recall that to < /°o 5 therefore we may replace po 
with to m (|5.39p . By (|5.38p and (|5.39p . we have that for every r, < r < hi, 

(5.40) \\ u \\ L2(n) <C(E + e)^J 

with C > and 7 G (0, 1) only depending on K, L, k, Mq, ^ and Setting 
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D \2 pj' 2D 1 * \E + e 

we obtain 

(5.41) ||«||r2 ( n) < C(E + E )(r* + t~°C), for every r, < r < r , 

with C > and 7 E (0, 1) only depending on K, L, k, M , & and ^. Denoting 

let us estimate from above inf <3?(r). To this aim, it is convenient to introduce a new 

0<t<t 

parameter I, related to r by 

/ 1 \ 



We compute 



*M = (i^r) + (i^r) «*>< 



log 7 



Let us temporarily assume I < 1, and let us use the inequality e s < 1/s with 

1 



C 

We obtain 



log 7, 



\ Z0 / \ l-Jfl+o-) 

1 \ / 1 



Us + Us 



Let us choose I = 1+ ^ +er , so that fo? = 1 — 1(1 + <r) = 1+ ^ +0 . , and < / < 1. Denoting 

M = min{W, 1 - 1(1 + a)} = ^ , 

1 + 17 + 0" 

we have 

$(t) < 2 



1 V 



x iogj; 

Now, if f < r , then min <J>(r) < 2 f r^r^ • Otherwise, if ( -r^—") > r , then 

_ o<t<t ~ V g c/ V g c/ 

min $(r) < 3>(to) < C, with C > only depending on K, L, k, M and 

0<r<ro ^° 

Furthermore, I r-^r I > T n\ an d hence min $(r) < 4 A -By (E3D and 

v log ?/ ~ v v log 7/ — 

recalling the definition of £, we arrive at 

/ 1 \" 

(5.42) <C(£; + £ ) 



"7 log IS, 
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where C > and fj,, < \i < 1, only depend on if, L, «, Mo, ^ and Therefore, with a 
possibly new choice of C, we obtain 



(5.43) 



u\\ L *(n) < C{E + e) 



log^ 



where C > and /i, < /i < 1, only depend on if, L, k, Mq, ^ and j^. 



□ 



6 Proofs of the stability results for the Cauchy problem 

In order to reduce the problem of continuation from Cauchy data to a problem of continu- 
ation from an open set, we shall introduce an augmented domain as follows. 

We recall that S is an open Lipschitz portion of dO, with constants po, Mq according to 
Definition 11.41 and it has size at least p\ > according to Definition 11.51 

We denote t] : [0, +oo) — > [0, 1] as follows 



(6.1) 



v(t) = { 



{ 1, 

4(i-t), \<t<\, 



t > i. 



lo, „^ 2 . 

Let PeS the point described in Definition 11.51 and let Z : B'p^ — » R be the Lipschitz 

Mo 

function appearing in (|1.16p for the local representation of O near P. According to (|1.16p 
we have 



Let us denote 
(6.2) 

Observe that 
(6.3) 



finu (p) = {x = {x',x n ) & r 



S^ 1 



x n > Z(x')}. 



Z-(x') = Z(x>) - ^ f ^T) , for every x> G B^. 

2 V Pi / M o 



\Z-{x')\ < \Z{x')\ + — , for every x' E P'_pq_, 

2 M 



(6.4) 

and therefore, since Mq > 1 



V 3; 'Z-(x / )| < \V x 'Z(x')\ + 2M , for every x' e B'^, 

M 



(6.5) ||Z" 

Next, we denote 
(6.6) 



+ Po||VZ- 



^4 = {x = (a/, x n ) e r po 



M 



^0 



I'll 



< P() H) + ( y + 2p ) M < 7 - P qMq. 



Z (x) < x n < Z(x ; )}, 
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(6.7) S = {x = (x',x n ) G r^ )P0 | \x'\ < ^-,x n = Z(x')}, 



(6.8) n = nus uA 

It is a straightforward matter to verify that 

i) if f2 has Lipschitz boundary with constants po, Mo, then has Lipschitz boundary with 

constants IMq. 

ii) using the coordinates employed in the construction we have 



(6.9) Q D T n pi ■ 

4M ' 4 

hi) if we denote, for any r > 0, 

(6.10) C7 = {x = (x, x n ) G M n | - r < x n < -M |x'|}, 

(6.11) C+ = {x = (x',x n ) G M n | M \x'\ <x n < r}, 
we also have 

(6.12) ADCJ,. 

4 

iv) if we denote 



pi + - 1 
(6.13) r = — 



Pi 



8M 2 



Y 7 1 + M5 2 + 1 



(6.14) Xo= ( ,r -^), 
we obtain that 

(6.15) B ro (x )cq 

8 

In the next lemma we continue to use the coordinate system centered in P and the 
notation described in Definition 11.11 and used above. 
We denote 



(6.16) T~ P1 n = {(x',x n ) G Tpi_ \x n < Z(x')} 

that is 



( 6 - 17 ) r ^ m = T ^, P1 \ a 
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Lemma 6.1 Let g G if 2 (£ n r_pi_ ). Then there exists v £ H 1 ^ P1 ) suc/i i/iat 

(6.18) ^|snr P1 =5 in the sense of traces, 

and 

where C > on/y depends on Mq and 

Proof This is a well-known fact. It suffices to prove it first in the reference situation when 
Z = 0, Mo = 1 and p\ = 1, see for instance [57, Lemma 6.9.1]. 

Next, by a scaling, and by our convention on norms as described in Remark 11.31 we 
obtain (|6,19p when Z = 0, Mq = 1 and p\ > is arbitrary. At this stage the constant C in 
()6.19p shall depend on ^ only. 

Finally, by a bilipschitz change of coordinates we may pass to the general (non-flat) 
case, at the price of admitting that C also depends on Mq. □ 

Let us now define 

{u, in Q, 

where u £ H 1 (f2) is a weak solution to the Cauchy Problem 11.61 and v is the function 
introduced in the previous lemma. 
We denote 

(6.21) fti = nu(snrp 1 „ )ur P1 . 

We obtain immediately that 

(6.22) uGH 1 ^) 
and also the following extension theorem. 

Theorem 6.2 (Extension) There exist f £ L 2 (Q,i), F £ L 2 (^i;IR n ) such that 

(6.23) ll71lL»(nx) + -wHinn^) < c e -±P- 

Po Po 

and u satisfies in the weak sense 

(6.24) div(,4V5) + cu = / + divF, in O x . 
Here C > only depends on Mq, K, k and 
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Proof Let ip be an arbitrary test function in Hq(Qi), with support compactly contained in 
fii. Evidently <p\ a e U E). 

Denoting for simplicity = , we compute 



M 



PI ' 



(6.25) - / (AVu ■ X7(p - cup) = - (AVu ■ V<p - cuip) - / (AVv ■ Vip - cv<p) = 
J six isi Jr~ 

= - ijtp+ / (f(p-F-V<p)- / (AVv ■ V<p - cv<p). 
Jt Jq Jr- 



Let us also denote Ei = Enr_Pi_ and let us set 



(6.26) *(<p) = [ = pr'^T I 1><P- 

Jt p Jt 

We have 

(6.27) 1*^)1 < pr'M^^ykA^^ < Cp^vMhIW 

Here C > is the constant for the trace imbedding Hq(Qi) <^-> i?5(Ei) which only depends 
on the Lipschitz character of Ei. Hence C only depends on Mq and 

Therefore * € an d its norm is bounded as follows 

(6-28) 11*1^-1^) < Cp n Q ~\ 

where C > is the same constant as above. 

By the well-known Riesz representation theorem in Hilbert spaces, we can find j\ G 
L 2 (fti), F x € L 2 (Oi;M n ) such that 

(6.29) Poll/ilU^no + ll^ill^cnuR") < -^11*11^-1(00 ^ 

Po ^° 

and 

(6.30) *(¥>)= / h<p-F!-V<p, for every y> e ^(fii). 

Note that the powers of po appearing above are calculated according to the fact that the 
appropriate scalar product for i7 1 (fii), as derived by our conventions on norms, is given by 

<¥>> < ^% 1 mi) = ~Ti f ff' + Po V V> • V^', 
Po ^Hi 

and, analogously, the L 2 (E)-scalar product is to be meant as follows 

(<P, <p')l 2 (e) = -^=r / </V- 
Po 

We define 



(6.31) / = 



f - fi, in ft, 
cv-fi, in r _ 
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(6.32) F 
We obtain 



F-F 1 , in ft, 
AVv-F u inT- 



(6-33) ll/l|jy«(ni) < II/IIl 2 (c) + ll/i||z^(ni) + -^IHl^t-), 

Po 



(6.34) H-^llL 2 (ni;M™) < H-F||L 2 (n ; R™) + 11-^1 llL 2 (ni;R") + -^11 II L 2 (r~ ;R") • 

Hence 

( 6 - 35 ) II/IIl 2 ^!) + — ll^llL 2 (ni;]R™) < 



< 



L 2 (Q) + — H-^llL 2 (n ; R")^ + ^II/iIIl 2 ^!) + — ll-Fl IIl 2 (Qi;R™)^ + \\v || H 1 (T~) ) 



where C > only depends on k and .fT. Thus, recalling Lemma [6. II and the assumptions on 
/, F, g and ip, we deduce (16.230 and from (16.250 we obtain (I6,24|) , □ 

Remark 6.3 It is evident how, by the above theorem, we can translate a Cauchy problem 
into a problem of propagation of smallness. The philosophy is as follows. Whenever we are 
given p\ > and P £ E such that Tpi_ (P) n E is a Lipschitz graph, and once we have 

shown that u\-p- = v is "small" in T~ Pl , then we can propagate such a "smallness" in a 

Holder fashion to any compact subset of T_pi_ . That is, we obtain that u = u\q is also 

"small" in a Holder fashion in a region near P inside ft. This rough argument is made 
precise in Theorem 1 1.7\ which we are going to prove here below. 

For technical reasons, we found it convenient to make use of the augmented domain 
ft C fti rather than fti itself. The advantage is that when ft is assumed to be globally 
Lipschitz, the global Lipschitz regularity of ft is proven in a more transparent way. This 
fact will be somewhat helpful later on for the proof of Theorem 11.91 



Proof of Theorem 11.71 Let us set xq and 7"o as in (|6.14p and in (|6,13p . respectively. Let 
h = min{2Co/90) if}? Co as in the thesis of Theorem 14.11 Notice that, by (|6.13p . h < 
Denote 

(6.36) G = GU T_£L ex- 

s,m ' 8 



Recalling (j6.9j) . it is easily verified that G C ft, it is connected and also 
(6.37) dist(G, 90,) > min { h, 1 = h. 



8Mn 



Moreover, by (|6.15p . 



(6.38) Bro(x ) cAnrjn. ex C G. 

2 8A/ ' 8 
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Hence we can apply Theorem 15.11 with u, f, F, O, G replaced with it, f, F, Q, G respectively. 
By Lemma 16.11 and by (j6.20j) we recall that 

(6-39) \\u\\L*(B ro (x )) < Ct], 

(6-40) \\u\\ L2{Vl) <C{E Q + ri), 

where C > only depends on Mq and Moreover, noticing that f2 D C^, it is immediate 
to estimate \Q\ < C|0|, with C > an absolute constant. Thus the thesis follows. □ 

Remark 6.4 Observe that in Theorem 1 1 . 71 the assumption u G H (ft) is not really neces- 
sary. This assumption is made just for the sake of simplicity. A more appropriate assumption 
would be 

(6.41) u£ p| H x (K). 

Kcnuz 

Remark 6.5 As is evident from the proof, the stability of Holder type could be obtained 
also on any connected subset G of f2 having a positive distance h from S' = dQ \ a (rather 
than from dfl). However, unless some additional assumption on the shape and regularity 
of is made, it might not be possible to specify in a precise manner how the constant C 
and the exponent 5 appearing in (|1.35p behave with respect to h as h — > 0. This is in fact 
an important issue in view of obtaining a global stability bound when no global regularity 
(Lipschitz) information on d£l is available. See Section [7] below for further discussion on 
this issue. 

Proof of Theorem 11.91 Similarly to what we did in the previous proof, we apply Theo- 
rem [5T3] with u, /, F, Q replaced with u, /, F, Q, respectively. As before, we consider the 
ball B ro (xo) as defined in ()6.13[) . (16. 14ft and we use the fact that, by Lemma IBTTI 

(6-42) \\u\\L2(B ro (x )) < Cr}, 

(6.43) \\u\\ m(ii) <C(E + r,), 

with C > only depending on Mq and □ 

7 A generalization of the stability results for the Cauchy 
problem 

Theorem 7.1 (Global stability for the Cauchy problem - generalization) Letu G 
H l {Q,) be a solution to the Cauchy Problem \l.b\ where S satisfies the conditions in Defini- 
tion\L]\and Definition\Tj& f G L 2 (M n ) and F G L 2 (M n ;R n ) satisfy (Tl~29l) and g G (S), 
tp G H~z (£) satisfy (jl.25p . Let us assume that there exists a family {G^}, < h < h, h 
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as in Theorem [777} of connected open sets Gh C f2 satisfying the conditions (jl.33p . (|1.34p 
and a/so 

(7.i) |n\G h | <Qp%(- 



/or given Q, $ > 0. //, given E > 0, p > 2, we a-priori assume that 
(7.2) IM|j>(0) < 



(7.3) 

where 



i L2ffi) < ( ^ +g+ ^( e( / + + ; +??) ), 



C 1 

(7.4) w(t) < jy, /or < f < -, 

W - (log|logt|) 5 ' 

where C > and S 1 , < 5 < 1, onZy depend on K , L, n, Mq, j^, Q, $ and p. 

Remark 7.2 In a completely analogous fashion, a global estimate of propagation of small- 
ness of log-log-type in a general connected open set O, could be stated and proved. 

Proof of Theorem 17.11 By Holder inequality, and by our convention on norms (Re- 
mark [L3]), we have 

i i 

(7-5) Il«lli3(n) < [j^J W u hp(n) < CE, 

where C > only depends on jp and p. Hence, by Theorem II. 7[ 

(7-6) \\u\\ LHGh) < C h (E + e + V ) [ E £ + + £ V +7 ^ ^ 

where 

i 

/lill \ 2 c 2 \n\ 
(7.7) C A = CifLiJ and 6 h > a T ( 

with a £ (0, 1) only depending on i^, L and k, Ci only depending on K, L, k, Mq and ^j, 
and C2 only depending on K. 
We notice that 



(7.8) \\ u \\ L2{Gh) < Ch (E + e + n)( E £ + + £ V + ^ 
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This last step has no substantial importance, it is made only because, at some later stage, 
it will be convenient to have that the ratio e ^^ +r ^ is strictly less than e _1 . 
Again using Holder inequality, we also have 

Setting 

(7 - 10) s= w 0<s - aa= w 

we obtain 

(7.11) \\u\\ L 2 (n) <C(E + e + 7])^h a ~ C2S +s°y for every sG (0, s ], 
where 



(7.13) ° = -(\- l -r 

y ' n \2 pj 1 

and C only depends on K, L, k, Mq, Q, -d and p. Note that also sq only depends 

Pi Pq 

on K, L, k, Mq, ^, A simple calculation (see for instance [5, Proof of Proposition 3.1]) 
gives 

(7.14) inf ( S -h a ' C2S +s D ) < C(log|logr|)- 5 , for every r G (0, e' 1 ), 
where C, S > only depend on K, L, k, Mq, ^, Q, $ and p. □ 

Pi Po 
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